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It is a pleasure for me to announce two new books: 

Our long time member Paul Drijvers is editor and contributor of 


SECONDARY ALGEBRA EDUCATION 

Revisiting Topics and Themes and Exploring the Unknown 

ISBN 978-94-6091-333-4 paperback 
SensePublishers, October 2010, 236 pages 


You can find a preview at: 

https://www.sensepublishers.com/product info.php?manufacturers id=68&products id= 
1 141&osCsid=27a01dc4ab7f32825afe0a4b9635b1c2 


I recommend visiting Paul's website: http://www.fi.uu.nl/-pauld/ 


And there is another one. DUG-Member Thomas Himmelbauer - he is author of several DNL- 
contributions - is not only an excellent mathematician, he also writes mystery novels. Just now 
his second book has appeared: 

TOD IM GYMNASIUM 

ISBN 978-3-902784-00-1, Taschenbuch, 205 Seiten 
Federfrei 2010 

His first book was TOD IN PANNONIEN (Death in Pannonia). The nice thing is that both stories 
are giving the atmosphere of the crime scene. TOD IN PANNONIEN is happening in the south- 
ern part of Burgenland where Thomas lives with his family. Much of the landscape and the peo- 
ple of this region can be found in his books. Southern part of Burgenland is “bilingual”, German 
and Croatian. 

TOD IM GYMNASIUM is laid in the atmosphere of an Austrian Secondary school. Reading it 
you can really smell the “Schulmief ’ (= “school stink”). 


Once again I recommend visiting Michael De Villiers ’ homepage. 

It was updated November 2010. 

http://mysite.mweb.co.za/residents/profmd/homepage4.html 

The link given below will lead you immediately to his latest newsletter containing many in- 
tersting links and lots of information. 

http://mysite.mweb.co.za/residents/profmd/newsletter.html 

One of the recommended sites will open a really enjoyable TILING SLIDE SHOW 

http://www.spsu.edu/math/tiling/tilings.html 
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bear DU(5 Members, 

Long, long overdue but now it is ready, 
bNL#80 can be downloaded. There are 
two reasons for the delay: The first one 
is our travel to Tansania in Novem- 
ber/bec-ember 2010 and the second one 
is a very extended and fruitful exchange 
of emails with Nils Hahnfeld in connection 
with the bEQME contribution. 



Lion in Serengeti 


Thanks for many Christmas and New 
Years wishes. Some of them showed nice 
graphics. They are included in this letter. 



Richard Schorn 


This bNL has only two - but very ex- 
tended - contributions. I didn't want to 
split them. 

We have an answer to a problem given in 
bNL#22 and a review of Nils' bifferential 
Equations tool for the TI-89, TT-92, and 
Voyage 200. 

Inspired by bEQME I tried to program a 
function for stepwise solving one type of 
bEs with bERIVE and TT-Nspire as well. 

I would like to put your attention to our 
Book-shelf (left page) and especially to 
the Modelling Books presented on page 4. 

All links recommended an pages 3 and 4 
have been checked and they should be 
valid. 



I received interesting letters from Roger 
Folsom and bietmar Oertel. They will be 
published next time. 



In July ACA 2011 will be held in Houston, 
TX. There will be an educational session 
and we will try to have again a special 
bERIVE and TT-CAS session. I will inform 
you as soon as possible. 


There is another conference in Germany: 
MNU in Mainz (7 - 11 April). Some bUG- 
members are giving lectures and work- 
shops there (R. Albers, K-H. Keunecke, 
W. Moldenhauer, P. Hof bauer, J. Bohm). 



www.bundeskonqress-2011.mnu.de 


Best regards as ever. 



Download all DTVL-Derive- and Tl-files from 

http : / / www . austromath . at/ dug/ 
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Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the DNL. It 
must be said, though, that non-English 
articles will be warmly welcomed nonethe- 
less. Your contributions will be edited but 
not assessed. By submitting articles the 
author gives his consent for reprinting it in 
the DNL. The more contributions you will 
send, the more lively and richer in contents 
the DERIVE & CAS-77 Newsletter will be. 


Next issue: March 2011 

Deadline 1 5 February 2011 

Preview: Contributions waiting to be published 

Some simulations of Random Experiments, J. Bohm, AUT, Lorenz Kopp, GER 

Wonderful World of Pedal Curves, J. Bohm 

Tools for 3D-Problems, P. Liike-Rosendahl, GER 

Financial Mathematics 4, M. R. Phillips 

Hill-Encription, J. Bohm 

Simulating a Graphing Calculator in DERIVE, J. Bohm 
Henon & Co, J. Bohm 

Do you know this? Cabri & CAS on PC and Handheld, W. Wegscheider, AUT 

An Interesting Problem with a Triangle, Steiner Point, P. Liike-Rosendahl, GER 

Overcoming Branch & Bound by Simulation, J. Bohm, AUT 

Diophantine Polynomials, D. E. McDougall, Canada 

Graphics World, Currency Change, P. Charland, CAN 

Cubics, Quartics - Interesting features, T. Koller & J. Bohm 

Logos of Companies as an Inspiration for Math Teaching 

Exciting Surfaces in the FAZ / Pierre Charland 's Graphics Gallery 

BooleanPlots.mth, P. Schofield, UK 

Old traditional examples for a CAS - what's new? J. Bohm, AUT 
Truth Tables on the Tl, M. R. Phillips 

Where oh Where is It? (GPS with CAS), C. & P. Leinbach, USA 
Embroidery Patterns, H. Ludwig, GER 

Mandelbrot and Newton with DERIVE, Roman Hasek, CZ & Rob Gough, UK 

Snail-shells, Piotr Trebisz, GER 

A Conics-Explorer, J. Bohm, AUT 

Tutorials for the NSpireCAS, G. Herweyers, BEL 

Some Projects with Students, R. Schroder, GER 

Dirac Algebra, Clifford Algebra, D. R. Lunsford, USA 

Treating Differential Equations (M. Beaudin, G. Piccard, Ch. Trottier) 

Structured Combinatorics, D. Oertel, GER 
Statistics with Tl-Nspire, G. Herweyers, BEL 
Cesar Multiplication, G. Schodl, AUT 

and others 

Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wiirmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag. Josef Bohm 


The DERIVE-NEWSLETTER is the Bulle- 
tin of the DERIVE & CAS- 77 User Group. 
It is published at least four times a year 
with a contents of 40 pages minimum. The 
goals of the DNL are to enable the ex- 
change of experiences made with DERIVE, 
77-CAS and other CAS as well to create a 
group to discuss the possibilities of new 
methodical and didactical manners in 
teaching mathematics. 

Editor: Mag. Josef Bdhm 
D'Lust 1, A-3042 Wiirmla 
Austria 

Phone: ++43-(0)660 4070480 

e-mail: nojo.boehm@pgv.at 
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The MacTutor History of Mathematics archive 

Biographies Index, History Topics Index, Additional Material Index, Famous Curve Index, 
Mathematicians of the Day 

http : // www- groups . dcs . st-and. ac . uk/ -hi story/ 


The German MathCad website 

http : //www . ptc . com/products /mathcad/ 


Wolfram Library Archive 

http: //library .wolfram. com/redir/ 


Eric’s Treasure Trove of Mathematics 

http : //mathworld. wolf ram. com/ 


Featuring over 2100 applications contributed by the Maplesoft user community 

http : //www . maple so ft . com/ applications/ 


Homepage of Cliff Pickover 

http : //sprott . physics . wise . edu /pickover /home . htm 


Topics in Mathematics 

In these pages, you will find links to various WWW resources on Mathematics. 
They are organized by topics. 

http : //archives . math .utk.edu/topics/ 


Mathematical Modules in Chemistry and Biology 

http : //science . kennesaw .edu/ ~mburke /modules/ 


The Geometry Center (University of Minnesota) 

Center for Computation and Vizualisation of Geometric Structures. The Geometry Center is now 
closed. This web site continues as a repository for much of the materials and projects from the 
Geometry Center. 

http : //www . geom .uiuc.edu/ 


Visual Index of all Uniform Polyhedra (R. E. Maeder) 

http : //www . mathconsult . ch/ showroom/ unipoly /unipoly 


The Math Forum - People Learning Math Together (Drexel University) 

http : //mathforum. org/ 

http : //mathforum. org/library/ 


Internet Center for Mathematics Problems 

http : //www . mathpropress . com/ 

with among others: 

http : //www . mathpropress . com/ archive /Rabinowitz Problems 1 9 63-2 0 05 . pdf 


Art from Code 
Enjoy the graphs 

http : //www . artf romcode . com/ 
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The Spanky Fractal Database 

http : //www . nahee . com/ spanky/ index . html 


Fractals - Chaos - Attractors a.o. 

http : / /local . wasp . uwa . edu . au/~pbourke/ fractals/ 


Yahoo-Science-Mathematics 

http : //dir . yahoo . com/ Science /mathematics 


Euclid’s Elements 

This edition of Euclid's Elements uses a Java applet called the Geometry Applet to illustrate the dia- 
grams. 

http : / /alephO . clarku . edu/~dj oyce/ j ava/ elements/ elements 


An interesting interview about “Darstellende Geometrie” as subject on Secondary schools (in 
German) can be found at 

http : //der standard . at/12 9557 0821131/Interview-Freie-Formen-f ordern- 
neue-geometrische-Modelle 


Download free e-books from 

http : / /bookbon . com/ uk/ student and http : / /bookbon . com/ de/ studium 

(New publication: Introductory Finite Difference Methods for PDEs) 


You are interested in Modelling? 

Then I can recommend Hartmut Bossel's 4 books: 

Systemzoo 1-3, Systeme, Dynamik, Simulation, Books on Demand, Norderstedt (German) 
System Zoo 1-3, Systems and Models, Books on Demand, Norderstedt (English) 

The books cover models of the following fields: Elementary Systems, Physics, Engineering, 
Climate, Ecosystems, Resources, Economy, Society and Development. They are based on 
the modeling software VENSIM PLE which can be downloaded free of charge. 

VENSIM PLE (free download for educational use) 

Ventana publishes Vensim which is used for constructing models of business, scientific, environ- 
mental, and social systems. 

http: //www. vensim. com/download . html 


This is nice tutorial for VENSIM PLE: 

System Dynamics Resource Page of the Arizona State University 

Among others you can find - and download a twenty-three page reference for Vensim PLE (pdf). 

http : / / www . public . asu . edu/~kirkwood/ sysdyn/ SDRes . htm 


There is also a CD available which contains all 100 models which are treated in the System 
Zoo books: 


Systemzoo, co:Tec: www . corec-verlag . de 
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It is very unusual but you will find the abstract as part of the article on page 9, Josef 


Using Rational Arithmetic to Develop a Proof 
“What Josef and Carl Saw” 

Josef Bohm, Wiirmla, Austria, and Carl Lewis Leinbach, Gardener, USA 
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Finding a Limit via 
Geometric Reasoning 

Carl Leinbach and Marvin Brubaker, USA 


Consider the following sequence of points: P n = < 


[0,0] 

n = 0 

[0,1] 

n = 1 

[1,0] 

n = 2 

V«-3 + 2 P n-2 

otherwise 


Notice that this sequence is defined recursively. DERIVE allows us to make recursive definitions. We 
use the IF statement. 


P (n) : =IF (n=0 r [0, 0],IF(n=l, [ 0 , 1 ] , IF (n=2 , [ 1 , 0 ] , 1/2P (n-3 ) +1/2P (n-2 ) ) ) ) 

In this case we had to nest the IF statements three deep. That is because we had three special cases. 
This function, because of its recursive nature, is slow to evaluate for an n of any size, whatsoever. 
Nonetheless, author 


VECTOR ( P (n) ,n, 0, 10) 

and plot the sequence. 
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The next figures show the evaluation of the first 10 terms of the sequence and also the first 20 terms. If 
we move the crosshair on the graph where the plot is dense, i.e., the point of apparent convergence we 
get a reading of approximately [0.4, 0.4]. 




We can zoom in and then we read off the coordi- 
nates of the crosshair [0.40029, 0.40042]. 

We can show the last term of the sequence given 
right above and we get a similar result: 
[0.40039..., 0.40039 ...]. 

Of course, we had not proved any result. How- 
ever, the visual evidence is convincing that a 
limit does exist ([0.4, 0.4]?) and we have a visual 
illustration of the process of convergence. 



FIRST(REVERSE (VECTOR (P (n) , 
FIRST(REVERSE (VECTOR (P (n) , 


0 , 

0 , 


20 ))) = 
20 ))) = 


' 205 205 ' 

. 512 ' 512 . 

[0.400390625 , 0.400390625] 


The challenge is still there: Proof that the limit is [0.4, 0.4]! 


Dear Carl, 


The History of the Lecture 

8 January 2010 


I am now revising DNL#22 which contains Carl's and Marvins' s Lab #2, 
"Finding a Limit via Geometric Reasoning". 

I had to change some things due to the fact that DERIVE has changed 
a lot since 1996. I attach the revised contribution. Hope that you 
are satisfied with the new form (including a small program) . 


My question is: do you have a proof for the limit [2/5, 2/5]? 

Best regards 
Josef 
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11 January 2010 

Josef - 

I have not started on Lab 2 , but hope to get to it before we 
leave on Wednesday morning. I have been working on meeting the (now 
revised) deadline before our Costa Rica trip. I enclosed te vastly 
revised paper in the hopes that you may find the example that I did 
on Time Since Death useful for your upcoming workshop. The referees 
wanted me to make my examples more "beefy", i.e. do some more sub- 
stantial mathematics and involve the CAS more than I did in the 
original paper we submitted. 

Dear Pat and Carl, 

please don't hurry - the proof is not so important. Enjoy your 
holidays . 

12 January 2010 

Josef - 

While I was in the doctor’s examination room waiting for the 
doctor to arrive, I tore off a piece of the paper covering the ex- 
amination table and started to write out terms of the sequence. I 
got up to 16 terms. 


5. Then prove that the lim(P(4*i)) = 2/5. 

At the moment everything is based on my suppositions, not proven 
fact. I will keep working. Just wanted to keep you up to date. 

some days later 

Dear Carl, 

Thanks for your efforts. 

I am on a very similar way - to investigate the pattern of the nu- 
merators . 

Hi Carl, 

I attach my ideas for proving the limit. 


27 January 2010 

Josef - 

I have attached the proof of the limit. I worked on it mainly 
on the plane ride to Costa Rica and a little bit during our visit to 
Costa Rica. It took a little more than I expected and as I note 
there is still one part that I want to clean up. I gave you an out- 
line of that part. It is essential to the argument and I don’t like 
the fact that it gets rather messy with the arithmetic. 
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3 February 2010 

Josef - 

I sent you this about a week ago and hadn’t heard back. I was 
wondering what you thought. I think that it could make a good talk 
on combining the use of the rational arithmetic display of DERIVE to 
stimulate conjectures for solid mathematical analysis and then de- 
veloping a proof. This is what we have been talking about for 
years. What do you think? BTW, I see that your though path and 
mine crossed at few crucial points. I was thinking that maybe we 
could develop a joint DNL article or a TIME talk on this type of use 
of DERIVE. Once again, what do you think? 

-Carl 


8 February 2010 

Josef I have mentioned a joint presentation at Malaga or a DNL arti- 
cle (your choice) . Here is how I thought it could go: 

History: The DNL #22 Article attributed to Marvin and Carl; a re- 
quest from Josef for an analytic proof of the limit 

Observation Phase : Writing a brief program to examine terms of the 
sequence; the advantage of the rational arithmetic calculations and 
print out of DERIVE (and other CAS's) 

Con j ectures : What Josef saw (even though we worked independently, 
you were first); what Carl saw; putting conjectures to the test: 
Using mathematical induction to construct a proof 

What do you think? I like the idea, because it uses a skill that we 
hope to develop amongst our students and uses CAS in much more than 
a "button pushing mode", which is what some of our antagonists ac- 
cuse proponents of using CAS in teaching say we are professing. 


16/18/20 February 2010 

Josef - 

Here is the promised draft of the Malaga presentation. Let me 
know what you think? Once we have the final form for the abstract, 
I will submit it. 

-Carl 

Dear Carl, 

It looks good, I am busy filling the gap(s) in my PROOF. Maybe 
that we could add one sentence about possible generalizations 
(changing the initial values, . . .) . 

I attach a DERIVE file containing a general form for creating 
our sequence of points together with a nonrecursive way to create 
the sequence with the requested lim. 

Josef 


D 
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8 March 2010 

To time20106ctima.uma.es 

Please, find attached in this mail a (lecture or workshop) proposal 
for the (ACDCA strand) (TI-Nspire and Derive strand) (Please, indi- 
cate the appropriate format and strand) . 

This is a Lecture Proposal for the TI-Nspire & Derive Strand 
Thank you, 

Carl Leinbach 

Abstract 

It all began with an article in DNL #22 entitled Finding a Limit via Geometric 
Reasoning authored by Marvin Brubaker and Carl Leinbach. In that paper the limit 
of a recursively defined sequence of points was found by connecting successive 
points with straight lines, thus creating a nested sequence of triangles that seem to 
converge to the point (0.4, 0.4). While editing an archival edition of DNL #22, Josef 
correctly pointed out that the paper did not really have a proof of the limit, only a 
collection of heuristic evidence gained by zooming in on the suspected limit. He 
wrote to Carl asking if he had a mathematical proof that sequence did, in fact, con- 
verge to its claimed limit. Both Josef and Carl began independent work on the prob- 
lem. Their initial step was the same. They each wrote a small DERIVE program to 
print out the first n terms of the sequence using the CAS’s rational arithmetic display 
of the points. After this their two approaches differed. 

In this presentation both Josef and Carl will discuss their approaches to constructing 
a proof that the sequence converges to its claimed limit, thus supporting the visual 
evidence. They will also discuss the value of using the Rational Arithmetic to sup- 
port the discovery of a strategy to accomplish their mathematical goal. If time per- 
mits, the presenters will investigate applying their approaches to other sequences of 
points. 


How Carl Attacked The Challenge 


Let’s suppose that a student had seen the Fibonacci sequence and the proof that the limit of the ratio of 
successive terms of that sequence converges to the “Golden Mean.” 


lim 

" F r. 


1 + ^5 
2 


This approach simply can not be mimicked. It leads nowhere. WHY? 

A next approach might be to try to visualize the terms of the sequence and look for some patterns. 
Suppose we try to familiarize ourselves with the nature of the sequence without using the features of a 
CAS, i.e. print out the decimal approximations to the sequence: 


p 10 
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001 0 0.5 0.5 0.25 0.5 0.375 0.375 0.4375 0.375 0.40625 0.40625 0.390625 0.40625 0.39S4375 0 . 39S4375 0.40234375 0.3984375 

0 1 0 0.5 0.5 0.25 0.5 0.375 0.375 0.4375 0.375 0.40625 0.40625 0.390625 0.40625 0.39S4375 0.3984375 0.40234375 0.3984375 0.400390625 


0,400390625 0,400390625 0.3994140625 0,400390625 

0,400390625 0,3994140625 0,400390625 0,3999023437 


What patterns do you see? 


Here Is What Carl Saw 


O^p 10 O.S^O.5 0.25 0.5 0.375 0.375 

0 /o Q.f^0.5 0.25 0 . 5^,375 0-375 0.4375 


0.< 

0 


0.400390625 0.400390625 0.399*140625 0.400390625 

/ / / 

0.400390625 0.3994140625 0.400390625 0.3999023437 


Observation 1: Every term of the first sequence lags one term behind the second sequence. Thus, we 
really only need to deal with one sequence. 


Proof: (Using the Principle of Mathematical induction) 

Base Case: Look at the terms of the sequence printed out above 
General Case: Assume the result holds for all k<n. Then 

Pfi,\ = ^(^-3,1 + Pn-2,l ) = ^(^z-4,2 + ^-2,1 ) 

Pyi-\,2 = ~^{Pn-4,2 + Pyi-1,2 ) = ^(^z-4,2 + ^-2,1 ) 


( 1 ) 

( 2 ) 


Where P Wj i designates the n - th term in the first sequence and V n 2 the same term in the second se- 
quence. The second equality in both (1) and (2) are a result of the induction hypothesis. 


[r 1 


0 0.5 0.5 0.25 0.5 0.375 0.375 0.4375 0.375 0.40625 0.40625 0.390625 

A -4 , 


0 1 0 0.5 0.5 0.25 0.5 0,375 0 


4375 0.375 0.40625 0.40625 0.390625 0.40625 


0.400390625 0.400390625 0.3994140625 0.400390625 

A 

0.400390625 0,3994140625 0.400390625 0.3999023437 


0.40625 

0.3984375 


0,3984375 

0.3984375 


0.3984375 

0.40234375 


0.40234375 

0.3984375 


0.3984375 

0.400390625 


Observation 2: P 4n ^ = P 4n ^ for all n = 0,1, 2, 3, . . . 

Proof: At the moment, it seems like the definition of the sequence is not going to get us to an obvious 
proof of this conjecture. 


Let’s see if something pops out by looking at the sequence in its rational number presentation. So 
let’s turn to DERIVE: 
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pts(n , m := 1/2 , pt) := 

Pnog 

pt := [0, 0; 0, 1; 1, 0] 
k := 4 
#1: Loop 

If k > n 
RETURN pt 

pt := APPEND(pt, [m-(pt4.(k - 3) + pt.j.(k - 2))]) 
k :+ 1 

#2: pts(20) 1 


Q 

1 

1 

1 

1 

3 

3 

7 

3 

13 

13 

25 

13 

51 

51 

103 

51 


2 

2 

4 

2 

S 

S 

16 

S 

32 

32 

64 

32 

128 

128 

256 

128 

1 

1 

1 

1 

3 

3 

7 

3 

13 

13 

25 

13 

51 

51 

103 

51 

205 

2 

2 

4 

2 

S 

S 

16 

S 

32 

32 

64 

32 

12S 

128 

256 

128 

512 


Observation 3: P 4i _ x 2 = P 4i 2 = P 4i+2 2 for all i = 1,2,3, . . . 
Proof: Assume that the result holds for all k < i 



by definition of the recursive sequence. 

The next to last equality was a result of the induction hypothesis. 


Finally, P 4i+2 2 = ~[P 4l2 2 + ^ 2 ) = P 4i 2 by the sequence definition and the first part of this proof. 


If we combine Observation 1 and Observation 3 we have the proof for Observation 2. Thus, the part 
of the “Geometric Reasoning” that states that the limit of the sequence of points lies on the line y = x is 
indeed correct. 


But: 


What is the value of the limit? 


Finding the Limit of {P 4it i } 


D U 1 
0 1 0 


0 

1 

1 

1 

1 

i 

3 

r 

J 

13 

13 

a 

13 

41 

ii 

m 

11 

204 

205 


204 


I 

2 

4 

2 

£ 

£ 

is 

£ 

JJ 

12 


32 

12-fl 

12£ 


1ZH- 

512 

512 

102* 

512 

1 

1 

1 

1 

J 

1 

7 

1 

11 

13 

n 


51 

■5-1 

1Q3 

£1 


HH 



EL5 

2 

2 

4 

2 

£ 

£ 


B 

32 

32 

bt 

32 

lit 

1ZB 


12B 

512 

512 

LC.’J- 

512 

ZW£ 


Finding this limit and then invoking observation 3 and one more observation, we can easily use a clas- 
sic £, S proof to show that the limit of this subsequence is the limit of the entire sequence. 

Looking at the sequence of first coordinates, we see that the even terms for i > 2 (remember I call the 
first term Pq ) have successive powers of two in the denominator . Here is a DERIVE program and its 
result to look at this sequence: 


p 
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Even_Terms(n r ev , pt) := 

Prog 

pt := [0, 0; D, 1; 1, 0] 
ev := [0 f 1] 
k := 4 
Loop 

If k > n 
RETURN ev 

pf := APPEND(pt, [1/2 ■ (ptj (k - 3) + pti(k - 2))]) 

If M0D(kf 2) = 1 

ev := APPENDCev, [ptik|l]) 
k :+ 1 

Even_TermsC35) 

[ 1 1 3 7 13 25 51 103 205 409 819 1639 3277 6553 13107 26215 

Of If — f f r f f f f p P f f f r f P 

L 2 4 8 16 32 64 128 256 512 1024 2048 4096 8192 16384 32768 65536 


Observation 4: P 2n = P 2 ,._ 1)1 + 


(-o' 




for all i = 1,2,3, 


and 


1 

2 


denotes the floor function. 


Proof: Once again, we will assume that the result holds for all k < i. 

We go back to the basic definition for the basic definition sequence, P m and work from there. 



This argument is laden with notation and not terribly instructive, so let’s give only an overview of how 
it goes: 


Break the attack into two cases: i even and i odd i.e. 2 i a multiple of 4 and not a multiple of 4. 
It is really the first case that we want, but need to prove it for all even terms. Basically, Ob- 
servations 1 and 3 get the P 2 / 3,1 term above to a previous multiple of 4 and then we work 
back up. The arithmetic gets messy and the exponents are a little hard to handle, but it eventu- 
ally all works out. Note that the sign change always takes place at the multiples of 4. As was 
mentioned: Observations 1 & 3 are the keys. 


Observation 5: 


Pm. 1 


= P 




, nr 1 

2-4‘ 


for i = 1, 2, 3, ... . ,and thus, 


4U h 2 - 4 ' 2 h 4 * 


Proof: This is just a matter of extracting the terms from Observation 4. 


Observation 6: lim P 4nl 

/?— »QO ’ 


= — and, thus, lim P, 


An, 2 


2 

5' 


Proof: We turn this one over to DERIVE: 

k 

1 « C-l) 2 

I 

2 k=0 k 

4 


5 
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Finally, we need only show that the sequences of first and second coordinates converge. We show 
that they are Cauchy Sequences of Real Numbers and use the fact that the Real Numbers are a com- 
plete metric space, i.e. all Cauchy Sequences converge. 


Observation 7: The sequences } and } are Cauchy Sequences. 

Proof: Let s > 0, Observations 1, 3, and 4 have shown that for any two adjacent terms in the interval 
from 4 i to 4(7+1) the absolute value of the differences are: 0,^-,^-,^-, respectively. Take 


the largest of these differences, and say that 


P 


# 5 


< 


1 £ 

Now, choose N such that for n> N , < — and 

2 h « 


r r -i 

# 5 


< — then if m, n > N we have 
4 


Ip -P = 

r», 1 m, 1 


P «.> P 4[ f ],l 


+ 


]•>) ( F "■* P 4[ f ].l 

Thus, P n ,\ is a Cauchy Sequence and hence converges to the same limit as P^\. 



2^ 


( 2^ 

P r , 


_ 

P r n 

4[f] 

’ 1_ + 


l W 5/ 


< £ 


The sequence P 4 , 2 is just one term ahead of P 4; , and, thus, also converges to 


Carl Is Finally Finished! 


How Josef Attacked The Challenge 

Mv first anni'oncli: 

This was the function for creating the visualisation, giving a sequence of points: 


PCrO 

If n = 0 
[ 0 , 0 ] 

If n = 1 
[ 0 , 1 ] 

If n = 2 
[ 1 . 0 ] 

1/2 -PCn - 3) + 1/2 -Pfn - 2) 
VECTORCPCn), n, 0, 20)’ 





1 

1 

1 

1 

3 

3 

0 0 

1 

0 

2 

2 

1 4 

2 

3 

3 



1 

1 

1 

1 

3 

3 

7 

0 1 

0 

2 

2 

4 

2 

S 

3 

16 

103 


51 

205 






25b 


123 

512 






51 


205 

205 






123 


512 

512 







7 

3 

13 

13 

25 

13 

51 

51 

16 

3 

32 

32 

64 

32 

123 

123 

3 

13 

13 

25 

13 

51 

51 

103 

8 

32 

32 

64 

32 

123 

123 

256 


I inspected the (sorted) numerators of the 1st components (the x- values) in order to find the 
pattern: 


p 
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Because of the recursive nature of the definition it needs a long calculation time finding the list of the 
first 100 numerators! The next function works iterative and is much faster: 

pts(n , m 1/2. pt) 

Prog 

pt :z [0 P 0; 0 r 1; 1 P 0] 

k := 4 
Loop 

If k > n 
RETURN pt 

pt := APPENDfpt, [m-(pti(k - 3} + pti (k - 2))]) 
k :+ 1 


SORT (VECTOR (NUMERATOR (k), k, (pts(60))ul)) 

[0, 0, 0, 1, 1, 1, 1, 1, 3, 3, 3, 7, 13, 13, 13, 25, 51, 51, 51, 103, 205, 205, 205, 409, 319, 319, 
319, 1639, 3277, 3277, 3277, 6553, 13107, 13107, 13107, 26215, 52429, 52429, 52429, 104357, 209715, 
209715, 209715, 419431, 333361, 333361, 333361, 1677721, 3355443, 3355443, 3355443, 6710337, 
13421773, 13421773, 13421773, 26343545, 53637091, 53637091, 53637091, 107374133] 

which is (for the first 60 fractions) without counting repeated appearances: 

[0, 1, 3, 7, 13, 25, 51, 103, 205, 409, 819, 1639, 3277, 6553, 13107, 26215, 52429, 104857, 209715, 
419431,838861, 1677721, ....] 


Starting with 7 we have always a package of 4 values containing the first and then three times the next 


value. I 
n= 4 

investigated the sequence of values from above starting with 
7 = 2 2 + 3 

n= 5 

13 = 

2-7-1 = 

2 3 +2-3-l 

n=6 

25 = 

2-13-1 = 

2 4 +2 2 -3-3 

n=l 

51 = 

2-25+1 = 

2 5 +2 3 -3— 2-3+1 

n = 8 

103 = 

2-51+1 = 

2 6 +2 4 -3-2 2 -3+3 

n= 9 

205 = 

2-103-1 = 

2 7 +2 5 -3-2 3 -3+2-3-l 

n= 10 

409 = 

2-205-1 = 

2 8 +2 6 -3-2 4 -3+2 2 -3-3 

w=ll 

819 = 

2-409+1 = 

2 9 +2 7 -3—2 5 -3+2 3 -3— 2-3+1 

n= 12 

1639 = 

2-819+1 = 

2 10 +2 8 -3-2 6 -3+2 4 -3-2 2 -3+3 


The elements of the sequence formed by the first row of P(n ) from above are the numerators divided 
by 2 n . 

I start with the elements with n = 4, 8, 12, . . . and try finding a general formula for the numerators. 


This was for me the real ’’funny part” of the problem! 
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# 1 : 


42 \ 


#3; 


# 4 : 


4-i - 7 i - I 4-k 2 i - ? 14 

2 f 3- I 2 - 3-2 - I 2 

k=0 k=0 

4-i + 1 

2 3 

5 S 

( 4-i 4 1 \ 

2 3 


VECTOR 


5 


4-i 

2 


S 


i, 0, 10 




7 103 1639 26215 419431 6710887 

16 256 1 4096 ' 65536 104B576 ’ 16777216 

1717986919 27187790695 439804651111 1 

4294967Z96 ' 68719476736 ' 1099511627776 I 


107374183 

268435456 


4 i + 1 

2 3 

— — — + — - 4-i 

#5: 5 5 3*2 2 


4*1 5 S 

2 



Derive simplifies expression #1 to a nice formula. Applying VECTOR I can check the correctness of 
expression #5 and in the last step the limit of the partial sequence is given - what can easily be calcu- 
lated without a CAS, of course. 


I repeat the procedure for elements with n = 5, 9, 13, ... and end again with the limit — . 

4.1-1 l - 1 4.k 3 i - 2 4-k 

#7: 2 + 3-2. X 2 - 3-2 . X 2 - 1 

k=0 k=0 


I can proceed in a similar way for the remaining elements of the sequence. 


For/? = 6. 10, 14, ... 


4-f 2 1 - 1 4-k l-l 4-k 

2 + 3-2 - X 2 - 3- X 2 

kzO kzO 


And finally for n = 3, 7, II. 15, ... 
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- 4-1 - 3 



2 2 




' 3 

51 

819 

13107 

209715 

3355443 

VECTOR 



i 0 10 

_ 







l 5 

5 

J 


. 8 ’ 

128 

2048 

32768 

524288 ’ 

8388608 


53687091 558993459 13743895347 219902325555 3518437208883 | 

13421 7728 ’ 2147483648 ’ 34359738368 ’ 549755813888 * 8796093022208 J 


r 

- 4-1 - 3 " 



2 

2 


2 

< 5 

5 J 


5 


2 

All partial sequences tend to the same limit, so the limit is 


Without doubt the CAS was a very strong support for my calculations and checking the results, but . . . 

I must admit that I was not really satisfied with my PROOF because I could not show that the pattern 
of the numerators and of the fractions as a whole will remain until infinity. 

Inspired by Carl's PROOF and by the fact that only natural numbers are involved I was quite 
sure that a proof by induction must be the right “recipe”. 

I used my formulae - which I had derived in the previous attempt -for generating a list of all fractions 
appearing in the sequence: 


#4: 


pUO 


-k 

2 


5 



-k - 1 

2 3-2 


5 5 



-k - 3 
3-2 


5 


2 


+ 


5 



1 

1 

3 

7 


2 

4 

8 

16 


13 

25 

51 

103 


32 

64 

128 

256 


205 

409 

819 

1639 


512 

1024 

2048 

4096 

VECTORCpCk), k, 1, 21, 4) = 

3277 

6553 

13107 

26215 


8192 

16384 

32768 

65536 


52429 

104857 

209715 

419431 


131072 

262144 

524288 

1048576 


838861 

1677721 

3355443 

6710887 


. 2097152 

4194304 

8388608 

16777216 


We have a table of the first 24 different fractions appearing in the sequence. 
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Then I started from the very beginning: 

I came back to the original sequences of the 1 st and 2 nd components. My consideration was that both 
components are created in the same way, then I could stick to only one of them and I chose the 
x-coordinate. Function pts(n) returns the first n first coordinates of the sequence of points. 

p ts ( 41 ) 1 


) 0 1 

0 

1 1 

1 

1 3 

3 7 

3 

13 

13 

25 13 

51 

1 — 1 
i_n 

103 51 




2 2 

4 

2 8 

8 16 

8 

32 

32 

64 32 

128 

128 

256 128 


) 1 0 

1 

1 1 

1 

3 3 

7 3 

13 

13 

25 

13 51 

51 

103 

51 205 


2 

2 4 

2 

8 8 

16 8 

32 

32 

64 

32 128 

128 

256 

128 512 


205 

205 

409 

205 

819 

819 

1639 

819 

3277 

3277 

6553 

3277 

13107 

13107 

512 

512 

1024 

512 

2048 

2048 

4096 

2048 

8192 

8192 

16384 

8192 

32768 

32768 

205 

409 

205 

819 

819 

1639 

819 

3277 

3277 

6553 

3277 

13107 

13107 

26215 

512 

1024 

512 

2048 

2048 

4096 

2048 

8192 

8192 

16384 

8192 

32768 

32768 

65536 

26215 

13107 

52429 

52429 

104857 

52429 

209715 ' 






65536 

32768 

131072 

131072 

262144 

131072 

524288 






13107 

52429 

52429 

104857 

52429 

209715 

209715 






32768 

131072 

131072 

262144 

131072 

524288 

524288 . 







I prepared another tool: I wanted to address each single element of the sequence, used the formulae 
p(k) from above and took in account the fact that it is better to consider packages of eight elements in 
a row instead of only four. 


el(n) := 

If n < 4 

[0, 0, 1, OLn 

If M0D(n, 8) = 5 v M0D(n , 8) = 6 v M0D(n, 8) = 0 
Cp(4- FL00R((n + 3)/8) - 3))|1 

#6: If M0D(n , 8) = 1 v M0D(n , 8) = 2 v M0D(n , 8) = 4 

Cp(4- FLOOR((n + 3)/8) - 3))i3 
If M0D(n , 8) = 7 

(p(4- FLOOR((n + 3)/8) - 3))i2 
(p(4-FL00R((n + 3)/8) - 3));4 


#7: 

VECTOR(elCk), 

k, 

41) 












#8: 


0 , 0 , 1 , 

o 

1 

1 

1 

1 3 

3 7 

3 

13 

13 25 

13 

51 


51 

103 



2 

2 

4 

2 ’ 8 ’ 

8 ’ 16 

8 

32 

32 ’ 64 

32 

128 

128 

256 



51 

205 


205 

409 

205 

819 

819 

1639 

819 

3277 

3277 


6553 

3277 



128 

512 


512 

1024 

512 

2048 

2048 

4096 

2048 

8192 

8192 


16384 

8192 



13107 

13107 

26215 

13107 

52429 

52429 

104857 

52429 

209715 






32768 

32768 

65536 

32768 

131072 

131072 

262144 ’ 131072 

524288 





Please compare with the first row of simplified pts(41) ' from above. 


Generalization of the problem 

For keeping the procedure more general I introduce the matrix ini which is the matrix defined by 
points #2 and #3; the first point is the origin and m = — by default ( m can be changed). 
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ptss(n, ini, pt , m := 1/2) 

Prog 

pt := APPEND([ [0 , 0]] , ini) 
k 4 

#9 : Loop 

If k > n 
RETURN pt 

pt := APPENDCpt , [m-CptjXk - 3) + pti(k - 2))]) 
k :+ 1 

As the first and second coordinates are following the same rule, it is sufficient to investigate only one 
of them. I am choosing the x- coordinates. 



r 

x y 




1 1 


ptss 

41, 



#10: 


x y 



k 

L 2 2 J 

/ 


# 11 : 


x x + x 2-x + x x + 2-x 4-x + 3-x 3-(x + x ) 6-x + 7-x 

121 2121 21 21 21 

0, x , x , , , , , , , 

1222 4 4 8 8 16 


7-x + 6-x 12- x + 13 • x 13 ■ (x + x ) 26-x + 25-x 25-x + 26-x 52-x + 51-x 

2121 21 212121 


16 32 ’ 32 64 64 128 


51-(x + x ) 102 -x + 103-x 103-x + 102-x 204-x + 205-x 205-(x + x ) 

21 21 21 21 21 


128 ’ 256 ’ 256 ’ 512 ’ 512 


410 -x + 409 -x 409 -x + 410 -x 820 -x + 819 -x 819- (x + x ) 1638 -x + 1639- x 

21 21 21 21 2 1 


1024 ’ 1024 ’ 2048 ’ 2048 ’ 4096 


For me it is important to double check the single steps of the procedure: 

Substituting [0,1] for x = [xi,x 2 ] results in the coefficients of x 2 which is the list of the 1st coordinates 
of the points: 


0 0 10 

1 

1 

1 

1 

3 

3 

7 

3 

13 

13 

25 

13 

51 

51 

103 


2 

2 

4 

2 

8 

8 

16 

8 

32 

32 

64 

32 

128 

128 

256 

51 205 


205 

409 


205 

819 


819 

1639 

819 

3277 

3277 

6553 

3277 

128 ’ 512 


512 

1024 


512 

2048 


2048 

4096 

2048 

8192 

8192 

16384 

3192 


13107 13107 26215 13107 52429 52429 104857 52429 209715 

32768 ’ 32768 ’ 65536 ’ 32768 ’ 131072 ’ 131072 ’ 262144 ’ 131072 ’ 524288 


I substitute [1,0] for x = [xi,x 2 ] for obtaining the coefficients of x\ (= 2nd coordinates of the points): 


#13: 


0 , 1 , 0 , 


205 


111 
2 ' 2 ' 4 

409 205 


13 3 

2 ' S ' 8 

819 819 


512 1024 512 2048 2048 

26215 13107 52429 52429 


65536 


32768 


131072 


131072 


7 3 13 

16 s' 32 

1639 819 


13 25 13 51 51 

32 64 32 128 ' 128 

3277 3277 6553 3277 


103 


51 


205 


256 128 512 

13107 13107 


4096 2048 8192 8192 16384 8192 32768 32768 

104857 52429 209715 209715 


262144 


131072 


524288 


524288 


#14: x2_p(n) := el(n) 

#15: xl_p(n) := el(n + 1) 
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I split the fractions into their summands try to proof the pattern of the coefficients by induction. 

Assume that the rule is valid until element xl_p(«) with mod(«,8) = 0; we would like to find element 
xl_p(«+l)by 1/2 *(xl_p(«+l-3)+xl_p(«+l-2)) = = l/2*(xl_p(«-2)+xl_p(«-l)). 

Then x l_p(«-2) with mod(//-2,8) = 6 and x 1 p(/?- 1 ) with mod(«-l,8) = 7 will - hopefully - give 
xl_p(«+l) with mod(«+l,8) = 1. 

First of all I perform a - successful - check again (for n = 40): 

209715 

#16: xl_p(40) = 


524288 


#17: xl_p(41) = (xl_p(3S) + xl_p(39)) 

2 


209715 


209715 


#18: 


524233 524288 

I copied function el (n) because am needing the subexpressions for the different cases of mod(w,8). 

el(n) := 

If n < 4 

[0, 0, 1, 0] j,n 

If M0D(n, 8) = 5 v M0D(n, 8) = 6 v M0D(n, 8) = 0 


#19: 


(p (4. FLOOR ((n + 


- 3))*1 


If M0D(n, 8) = 1 v M0D(n, 8) = 2 v M0D(n, 8) = 4 
(p(4 ■ FL00R((n + 3)/S) - 8))|8 
If M0D(n, 8) = 7 

(p(4 ■ FL00R((n + 3)/S) - 3))j,2 
(p(4- FLOOR ((n + B)/S) - 3))j,4 


Then xl_p(^-2) (with mod(^,8) = 6): 

n + 3 


SUBST 


# 20 : 


# 21 : 


4- FLOOR 


8 


3 | J , n , n - 2 
1 


3 - 4 ■ FL00R(n/8 + 1/S) 

2-2 2 


#22: VECTOR 


3 - 4 ■ FLQ0R(n/S + 1/8) 

2-2 2 

+ , n, 8, 48, 8 

3 3 


1 13 205 3277 52429 S3SS61 

2 ' 32 ' 512 ' 8192 ' 131072 ' 2097152 


Expression #22 are elements #6, 14, 22, 30, ... 
xljp(n-l) (with mod(^,8) = 7): 


#23: 


#24: 


SUBST 


ff r 

" n + 3 ^ 

Y| ^ 

p 4- FLOOR 


- 3 n n - 1 

. 8 j 


L L 



2 j 


2 - 4 ■ FL00R(n/8 + 1/4) 

2 -3-2 


#25 : VECTOR 


f 2 - 4 ■ FL00R(n/ 8 + 1/4) ^ 

2 2 -3-2 

n fi /|0 fi 


" 1 

25 

409 

6553 

104857 

1677721 ' 

.5 5 j 


. 4 

64 

1024 

16384 

262144 

4194304 . 


Expression #25 are elements #7, 15, 23, 31, ... 
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This is - should be - the next element in the sequence xl_p(?z+l) (with mod(7z,8) = 1): 


SUBST 


#26: 


#27: 


n + 3 

p| 4 ■ FLOORl | - 5 | | , n, n + 1 


1 - 4 ■ FLOOR (n/S + 1/2) 

2 2-2 


#2£: VECTOR 


r 

2 

1 - 4 ■ FL00R(n/£ + 1/2) 

2 -2 

> 

n, £, 4£, £ 

j 


" 3 

51 

£19 

13107 

209715 

3355443 " 

^ 5 

5 


. £ 

12£ 

204£ 

3276£ 

5242££ 

£3££60£ . 


The next check holds: 


#29: 


#30: 



" 1 

13 

205 

3277 

52429 

£3 ££61 " 


" 1 

25 

409 

6553 

104 £57 

1677721 T 

2 l 

. 2 

32 

512 

£192 

131072 

2097152 . 


. 4 

64 

1024 

163£4 

262144 

4194304 _L 


51 £19 13107 209715 3355443 


L £ 12£ 204£ 3276£ 5242££ £3££60£ J 

Now follows the interesting step: 1/2*(#21 + #24) = #27 ?? 


# 51 : 


#52: 


1 

r 3 - 4 ■ FLOOR (n/S + 1/S) 

2-2 2 

r 

2 

2 - 4 ■ FLOOR (n/S + 1/4) V 

2 .3.2 

2 

h 

l ^ 

h 

l n 
J 

^ 5 

LTl 


2 - 4 - FL00R(n/S + 1/S) 

2 -2 


1 - 4 - FL00R(n/S + 1/4) 

2 -5-2 2 

+ — 

5 5 


DERIVE does not simplify further because it has no information about the nature of n. But we have: n 
is divisible by 8 (mod(7z,8) = 0). 


We know that: for all n with mod(7z,8) = 0: floor(7z/8 + 1/8) = floor(7z/8 + 1/4) = 72/8, so we can pro- 
ceed: 


2 _ 4- (n/S) 1 - 4- (n/S) 

2-2 2 -5-2 2 

#55: + 

5 5 5 


(2 - n)/2 

2 2 

#54: 

5 5 


We can do this with the CAS, too. Let's take in account that n is divisible by 8. I substitute n by 8 n 
and try again simplifying the expression. 
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1 

r 3 - 
2 -2 

4 ■ FLOOR (n/S + 1/S) 

2 


r 

2 


2 - 4 ■ FLOOR (n/S + 1/4) 

2 -3-2 

v 

2 


5 

5 


. 5 


5 ^ 



r 3 - 

4-FL00R(S.n_/S + 1/S) 




r 

2 - 4.FL00R(S.n_/S + 1/4) Y| 

1 

2 -2 


i 

2 

i 

2 

2 -3-2 


2 


5 

t 

5 

T 

. 5 

5 

Jj 


2 - 4 ■ FLOOR (n_ + 1/S) 1 - 4.FL00R(n_ + 1/4) 

2-2 2 -3-2 2 


5 5 5 

n_ :e Integer (0, ®) 

1 - 4-n_ 

2 2 


5 5 


1 - 4 ■ (n/S) 

2 2 


5 5 


(2 - n)/2 


I “simplify” expression #27 in the same way: 


1 - 4- FLOOR (n/S + 1/2] 
2 2-2 

#27: 

5 5 


1 - 4 ■ FL00R(n/S + 1/2) 

2 2-2 

# 35 : 

5 5 

1 - 4 ■ (n/S) 

2 2-2 

#36: 

5 5 

(2 - n)/2 

2 2 

#37: 

5 5 



We can repeat the procedure for all cases and proof show the identities of 

xl_p(^7+l) = l/2-(xl_p(^-2)+xl_p(^-l)) for all positions of n within a package of 8 in a row. 

It is obvious that for the second part of the x- value = x2_p(?z) the proof will also hold. 


2 2 

What we also can see is the fact that the full x- value will be — x x + — x 2 + f\{ri)x x + fl(n)x 2 where fl 
and fl are functions with T in the denominator. The same is happening with they- values. 
Calculating the limits, the functions are tending to 0 and the limit of the sequence of points with 

r 2 2 

[x 0 ,yo] = [0,0] will end in - (x, + x 2 ), - (y, + y 2 ) 


See an example: Initial points are [0,0], [5,-4] and [1 1,9]. 
If my idea holds then the sequence should end in 


|(5 + ll),|(-4 + 9) =[6.4,2]. 
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#47: 


#4S: 


r 

f 

i 

l n 

>> 

ptss 

100, 

_ 11 9 _ 

Jj 


6.4 2 

6.4 2 

6.4 2 

6.4 2 

6.4 2 

6.4 2 


[95, 


100 ] 


I introduce a more general function including variable pt (= a matrix) for the initial points: 

ptsG(n, pt, m := 1/2) := 

Prog 
k := 4 
Loop 

#49: If k > n 

RETURN pt 

pt := APPEND(pt, [m-CpttCk - 3) + pt j, fk - 2))]) 
k :+ 1 


#50: 


#51: 


#52: 


r 

r 

" 0 

0 

ptsG 

100, 

5 

-4 



_ 11 

9 

" 3602S79701S963S9 


L 562949953421312 
[6.4, 2] 




1JJ 

100 

1407374SS355327 " 
7036S744177664 _ 


Initial points are [-3, 5], [5, -4] and [11, 9]. What is the convergence point now, if there is one? 


#53: 


#54: 


ptsG 


100 , 


-3 5 

5 -4 

11 9 


100 


20406935 S115 225 16SSS49S60263931 


351S43720SSS32 562949953421312 


#55: 


[5.S, 3] 


Can you find out the rule? 


0.4-x + 0.4-x + 0.2-x , 0 .4-y + 0.4-y + G.2y 1 

3 2 1 3 2 lj 


Proof this! 


What happens if m * 1/2? Conjectures? Proofs? 
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So What Was The Role Of The CAS? 


© Although the DNL #22 Article said that it was finding a limit, it really only 
gave our intuition a “nudge”. 

© To really know that (0.4, 0.4) is the limit, a proof was required. The CAS can 
not construct a proof. There is no button to push. 

© This is where a “partnership” develops. The student, and instructor, have to 
understand what it is that the CAS and other technologies can do to help with 
the reasoning process. 

0 Visualization is a powerful aid. Sometimes it takes the form of graphical 
displays, other times it may be just to generate a large number of terms or 
examples, or, as it this case it was to give a display that made certain patterns 
“stick out.” 

@ As instructors, we need to “let a thousand flowers bloom”, i.e. let our students 
try their own strategies and exercise the limits of the CAS and other techno- 
logies. Our role is to gently critique and offer guidance through suggestions. 
In this case, a real strategy did not emerge until it became clear that the deno- 
minators were powers of 2. Everything else emerged from this very simple 
observation. 


Postlude 

Here’s What Rudiger Saw 

The original article was from DERIVE Newsletter #22. Rudiger Baumann sent a short note for the 
DERIVE Newsletter pointing to the fact that little generalization leads to Edward Sawada’s “Mis- 
guided Missile” contribution (also from DNL#22). 

Rudiger recommended the ITERATES-procedure because the recursive procedure is too slow. 

#35: pts_baum(r, 5 , ini, n) := ITER.ATES([b, c, r-a + s-b], [a, b, c], ini, n) 

This is the "Leinbach-Brubaker Sequence": 



f 

0 

0 


#36: pts_baurn 

0.5, 0.5, 

0 

1 

■ 50 



_ 1 

0 _ 

J 


And this is Edward's missile: 


0.9 , 0.1, 



#37 : pts_baurn 
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Playing with the parameters in Rudiger’s function leads to interesting patterns (limits?) 



f 

"00" 

\ 

#33: pts_baum 

0.15, 0.9, 

0 2 

, 50 



. 20 . 

J 


f 

"00" 

'N 

#90: pts_baum 

0.05, 0.9, 

0 2 

, 50 



. 20 . 

J 


#S9 red and #90 blue 


Let's produce a TWIN 



f 

"03" 

> 


f 

"30" 

\ 

pts_baum 

0.02, 0.3, 

1 1 

, 30 

, pts_baum 

0.02, 0.3, 

0 3 

, 30 



_ 3 0 _ 

j 



_ 1 1 _ 

J 




Again working with my “beloved” sliders: 

Why not trying to introduce sliders for the parameters r and s and investigate their influence on the 
sequence of points? 


Are this really conics? 


Try a proof!! 


If you find other (better?) proofs for the pre- 
sented problems then please send them. 


• 2.00 





IS] 

°00 

o 00 




Carl and Josef 
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Differential Equations Made Easy (2) 

Review and Description of some features, Josef Bohm 

In DNL#74 I presented Nils Hahnfeld’s tool DEQME 
and was busy with Menu FI dealing with 1 st order 
ODEs. In this article I’d like to proceed to Menu F2 
which offers treating 2 nd order differential equations. 

I’ll try to address all Menu options from 1 : through B: using the occasion to demonstrate parallel how 
to apply DERIVE for solving the problems and comparing with the solutions given by other CASs. 

I don't hesitate to admit that I learned - again after many years - about specific DEs - and I liked to 
do this. Most of the examples are from a textbook Differential Equations [1] . 

I'll start with 1 : Any 2. Order DE and observe the reaction of DEQME. My first package of differen- 
tial equations is: 

(1) y" = x 2 +3y- 6; y( 0) = 4,/(0) = -4 

(2) y" + 4y = l2x; y(0) = 5, y' (0) = 7 

(3) y' " + y = 2 sin(x) • sin(2x) 



Example (1) 




In order to make comparing the tools easier I will do that example for example. I will am starting with 
DERIVE and proceed with WIRIS, wx Maxima, MuPAD , and with TI-Nspire , too, of course. 

DERIVED Online Help informs about the syntax for solving 2 nd order ODEs: 

DSOLVE2(p, q, r, x, cl, c2) simplifies to an explicit general solution of the linear second 
order ordinary differential equation 

y" + p(x)-y' + q(x)-y = r(x) 

DSOLVE2_BV(p, q, r, x, xO, yO, x2, y2) is similar to DSOLVE2, but simplifies to a specific 
solution that satisfies the boundary conditions y=y0 at x=x0 and y=y2 at x=x2. 
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DSOLVE2_IV(p, q, r, x, xO, yO, vO) is similar toDSOLVE2_BV, but simplifies to a specific 
solution that satisfies the initial conditions y=yO and y-vO at x=xO. 

2 

#1: DS0LVE2_IV(0, -3, 6 - x , x, 0, 4, -4) 


# 2 : 


■/3-x. 


26 2--/3 


■J3-Y 


26 


3-x - 16 


v. 9 3 J 

Now let WIRIS try the job: 

probl : = solve(y ,, (x) = x ? +3y(x)-6,y(0)=4,y , (0) = -4);| 
probl 




I applied a “trick” introducing probl in order to avoid presenting the equation together with the result 
in one line which would have been difficult to print it here in a reasonable size. 

Next in the row is wxMaxima : 

(Sil) eqn_l: T diff(y,x,2) = x A 2+3*y-£; 

(Sol) 7=3 /+k -6 

a, 2 


(Si2) ode2(Sol, y, x) ; 


-J7V 3x " 16 

(So2) +%k2%e ^ 


(Sil7) ic2(So2, y=4 , T dif f (y, x) =-4 ) ; 

I 1 -\3 x } 1 --\3 x 2 

( 6-^3 " 10 ) "se * (6-^/3 + 10 ) %e * 3 x -16 

(Sol7) Y= + 


What about MuPADl 


ivp: =ode ( {y T T (x) =x A 2+3*y (x) -6, y (0) =4, y T ( 0 ) =- 4 } , y (x) ) 
ode( ^(0) = 4. y(0)= -4, -3-y(x)-x 2 +£j-y(x) + 6'?,y(x) 


solve (ivp) 
W3- 


2-V^ 10 

3 9 


VT 


2^/3_ 10 
3 9 


zl+Ii 


Last but not least see the TI-Nspire : 


deSolve(y" ] =x^+3y-6 and>{o)=4 andy(o)=-4 J x Jt y J 


LjT 1Q ) e E' x )'e'''' 3 ' X x 2 


y= I — - — h |-e 

3 9 


16 


3 9 


TI-Nspire works like the TI-92PLUS and Tl-Voyage 200. Please wait, DEQME is able to do a lot 
more than only solving this kind of DEs assisted by a nice form for entering the data. 


All systems behave the same so far. I can assure that things will change! 


D 
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Example (2) 


Fi 

r Ent4K ony 2. DKd4K DE. '' 

1 


| 


DE: |y 1 1 +4y=12x 1 


Exl: y 1 l= x*y 

Indep. uar.=: lx 

Dep. uar.=: |y 

OPTIONAL: Initial condition 

at x or t or . . =: |0 

y or s or . . =: 15 

y 1 or s 1 or . . = : 17 

^CEnter=0K !> CESC=CflNCEL^ 

MAIN RAD AUTD FUNC i/30 



General solution: 
y = 057 ■ cost 2 ■ x) + 058 ■ sin(2 ■ x) + 3 ■ x 
Using Initial condition, 
particular solution is: 
y = 5 ■ cos(!2 ■ x} + 2 ■ sin(2 ■ x} + 3 ■ x 


MAIN 


RAD liUTD FUNC 1/50 IflIH.I 


DERIVE: 

DS0LVE2_IV(0, 4, 12-x, x, 0, 


2 

B, 7) = 4-SIN(x)-C0S(x) - lO-SIN(x) + 3-x + B 


WIRIS: 

prob2 : = solve(y ,, (x)+4y(x) = 12x,y(0)=5,y , (0) =7) ; 

prob2 — Hy(x) =6-x-sin(2-x) ? + Sin ^ +6-x-cos(2-x) 2 +5-cos(2-x) 

sin (2 * x 1 

prob2 : = simplify(S-x-sin(2-x) 2 + y - — -+S-xcos(2x) 3 +5-cos(2-x) ); 

sin f2 ■ xi 

prob2 — 5cos(2x) + — -+6x 

WxMaxima : 


(*i20) 

(%o20) 


eqn_2 : T dif f (y, x, 2 ) +4 *y=12 + x; 

,2 


-y + 4 jr = 12 x 


d x 


(Si21) ode2 ( % o2Q , y, x) ; 

( % o21) y = %kl sin( 2 x ) + %k2 cos(2x)+3x 


(*i24) ic2 (*o21, x=0, y=5 , T dif f (y, x) =7 ) ; 
( % o24 ) y = 2 sin( 2 x)+5 cos(2 x ) + 3 x 


MuPAD provides the same result and so do the TIs! 


WIRIS needs an extra simplification - it does not simplify 6x sin 2 (2x) + 6x cos 2 (2x) to 6x. The WIRIS 
result seems to be wrong! I am using WIRIS itself to check the result - and the right side gives 24x 
instead of 12x. 


y"(x)+4-y(x) — *■ 24 -x 


1 ' 2 


I could have used also DEQME ' s 2nd order DE Checker to check the presented (correct) solutions! 
Example (3) (including an additional IVP): 



r Ent*r onv 2 . Drdtr DE. 

1 

DE: |y 1 1 +y= 2 *sin{x)*sin{ 2 x ; 1 1 

1 

Exl : y 1 1 =x*y 

Indep. uar.=: lx II 

Dep. uar.=: |y I 

OPTIONAL: Initial condition 

at x or t or . . =: I 

y or s or . . =: 1 

y 1 or s 1 or . . =: 1 

, CEnter= 0 K 3 CESC=C ANGELA ^ 

TVFE + [ENTER]=DK AND [ESC]=CANCEL 
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ifl^Y FF? YTF^t fht y fl t Fe^ t "i 
^—|fllgebra |Ca leather- |Frgro IQ |c lean Up| [ 


4 (@51 - (sin(x)) 4 ] ■ cos(x) + + G52j ■ sin(x) 


re si 


Pi m,.i 1 

General solution: 

-sin(x) ■ sin(4 ■ x) ,f ncc , . ^ v -,4l l 
y= g + [G55 - (sin£x)) J-cos* 

Using Initial condition, 

particular solution is: 

-sin(x) ■ sin(4 ■ x) ,f n / . v v -,4l l 

y = g + [2 - (sin(x)) J ■ cos(xr 


FUNC 1/50 


UTTET 


Working with DERIVE is interesting: 


EnUr an!- 1 £. Dk-J-jk DE. 


DE: |y 1 1 +y=2*sin<x>*sinC2x>~ 


n 


Exl : y 1 1 =x*y 
Indep. uar.=: [x^ 

Dep. uar.=: 


OPTIOHflL: Initial condition 
at x or t or ..=: Ijix2 

y or s or ..=: II 


y 1 or s 1 or . . =: I ~1 


MAIN 

RAD AUTD 

____ — -■ 

FUNC 1/30 


Pi=t 

sura 

*n 


General solution: 

y = -sinOO^sin^x) + ( e55 _ ^^4) . ^ 

Using Initial condition, 
particular solution is: 

^ +(2 -(sin(x» 4 ) cos(x) +(^- - 71 4 4 j sinCx) 

MlilN Fill [i liUTD FUHl 1/50 


DSQLVE2_IV 


0, 1, 2- SIN(x) ■ SIN(2-x) , x, , 1, -1 

2 


C05(x) -C0S(4-x) SIN(x) ■ SIN(4-x) COS(x) -C0S(2-x) 13-C0S(x) 

- + + + 


x 

f 2 


+ 1 


■ SIN(x) 


1(4 -x + 1) ■ C0S(ATAN (2 ■ x) - x) 


2 


13 


n - 32 -tt + 233) ■ SIN 

ATAN 


- X 



^ 2-(jt - 4) j 

J 


C05(3-x) 


S 


[Trigonometry := Expand, Tn'gpower := Sines] 
2 


C05(x)- 


2 - 


SIN(x) 


+ 1 


SIN(x) 


There are three different appearances of the result depending on the Trig Mode Settings. You need 
some skills (or intuiton and luck - just try) to find the appropriate settings in order to obtain a compact 
form of the solution. Plotting all solutions on the same axes and checking by substituting the solution 
into the given equation shows the identity of the expressions. 

SIN(x) ■ SIN(4-x) 


y(x) 


s 


+ (2 - SIN(x) )-C0S(x) + 


TT — 4 


k 2 


SIN(x) 


y rr (x) + y (x) - 2-SIN(x)-SIN(2-x) = 0 

WIRIS presents its solution in a very “extended” form. Plotting and applying the appropriate Trig 
Mode in DERIVE confirms the identity with the other results. 

solve(y , '(x)+y(x)=2sin(x) sin(2x),y(tt/2)=1 ,y'(7i/2) = -1) 


n.. r _,_ sin(x) -sin(-2x) 

sin(x) -sin(4-x) 

sin(x) -sin(2-x) ( x n \ 

Ll y(X) " 4 

8 

4 + l 2 ” 4 +1 j 


y(x) 


sin(x) sin(-2 x) sin(x) sin(4 x) sin(x) sin(2 x) f x n \ 

: = — 1 — 1 — 1 -+ + 1 sinfx) - 

4 8 4 1 2 4 J 1 } 


cos(4 x) cos(x) cos(x) cos(2 x) 13- cos (x) 


8 2 8 

This was the first step. Simplifying this result once more delivers the result as above (DERIVE): 
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simplify! 


_ sin(x) -sin(-2-x) _ sin(x) -sin(4-x) _ simi 
4 8 




WxMaxima (correct result): 


cos (3 x) + x 3±n[ x) - cos ( x) (%pi-4)sin(x) 

y = + 2 cos ( x ) 

S 4 


I believe that you are not surprised now that MuPad also delivers another form of the result - which 
proved to be wrong (satisfies the IVs, but does not satisfy the DE). 


solve (ivp3) 
23 ■ cos(x) 


16 


. , , cos(3 ■ x) cos(5 ■ x) 
+ sm(x) + ^ — -+■ ' ' 


32 


32 


71 ■ sin(x) + x ■ sin(x) 


Finally I apply tCollect on the DEQME - output 
and receive a new appearance (right). 

Then I substitute into the given equation and 
calculate the difference of the left side of the 
equation and the expected right side and hope 
that the result would be 0. 


fEiTmOY fit t f::t 'i 

| t |fl 1 gebra |Ca 1 c 

f FHt 1 

Other 

PrgmloTciean UpT I 

- sin(x) ■ sinH 

i-x) r 

2 - ( s i n( x)) ^) ■ cosifr 

y " S 

+ [ 

■ right(res) 

" s i n( x) ■ s i n( 4 ■ x) 

-+{ 2 - 

(sin(x))^] ■ cos(x) 

S 

■ tCol lect( rights res)) 
cost; 3 ■ x) + 15 ■ cos(x) + 2- (2-x-Ji + 4)sin( k 

8 f 

.„cos<x>+2#<2#x-n+4}#sin<x>>/8 

MAIN RAD AUTD 


FUNC 15/30 


■Fi-TmnY fT^ Y f:=t Y fht y ft Y Fir Y 'I 
t f — |fl 1 gebra |Ca 1 c |dther |Pr gn 1 0 |c 1 ean Up | [ 



(rights res)) + right(res) - 2 ■ sin£x) ■ 


"cos(x) ■ cos(4 ■ x) + 2 ■ sin(x) ■ sin(4 ■ x) - 2 ■ 

■ tCol lect( "cos(x) ■ cos(4 ■ x) + 2 ■ sin(x) ■ sin^ 
0 

t Collect Cans 

MAIN RAD AUTU FUNC l?,' 3* 


Both next problems read originally as follows: 

Verify thatjq andy 2 are solutions of the DE. Then find a particular solution of the form 
y = Ci yi + c 2 y 2 that satisfies the given initial conditions. 

(4) *y - 2*/ + 2y = 0; y x = x, y 2 = x 2 ;y( 1) = 3, /( 1) = 1 

(5) x 2 y" + 2xy'-6y = 0; y x = x 2 , y 2 = x" 3 ; y(2) = 10, y'{2) = 15 

Example (4) 


I try solving the DE and leave the check of the solutions y i and r-> for later. 


a 

^ EnUr qtiv £. Urdtr DE. 

P 

DE: |y 1 1 *x^2-2x*y 1 +2y=0 1 

Exl : y 1 1 =x*y 

Indep. uar.=: lx II 

Dep. oar.=: |y II 

OPTIONAL: Initial condition 

at x or t or . . =: 11 

y or s or . . =: 13 

y 1 or s 1 or . . =: 111 

, tEnter=0K !) £ESC=C ANGELA, 

MAIN RAD AUTD FUNC 1/30 
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Again we don’t encounter any problem using DEQMEl 

The next page shows the DE solved by DERIVE , WIRIS, wxMaxima , and MuPad\ 


D50LVE2_IV 


2 

x 


WIRIS : 


2 

. 0. x, 1. 

2 

x 


1 


= B-x - 


2 

2-x 


J 


Example (4) 

solvefx- ■y ,, (x)-2x-y'(x)+2y(x)=0,y(1 )=3,y'(1 ) = 1 ) — {{y(x) = -2 e 2 ln <*> +5 x}} 
j simplify(y(x) = -2-e 2 ' ln < x > +5-x) — * y(x) = -2-x 2 +5-x 

WIRIS needs again an extra simplification for e lnx = x. Then see wxMaxima . . . 


eqn_4: T x A 2* T diff (y f x, 2) -2*x* ’ diff (y f x, l)+2*y=0; 


— ry 


2 x 


— y 

d x j 


+ 2 y = 0 


ode2 ( %o ll f y f x) ; 
y = %kl x 2 + %k2 x 


ic2 ( %o 12 f K=l f y=3 f Miff (y f x) =1) ; 

2 

y = 5 x - 2 x 

. . . followed by MuPAD (compact form of the procedure performed on page 26): 
solve (ode ( { x 2*y T T (x) -2*x*y T (x) + 2 *y ( x ) = 0 , y(l)=3 , y T (!)=!} ,y(x) ) ) 
{5 -x — 2 - jc 2 } 

This runs pretty well, let's try the next example which looks very similar. 

Example (5) 



p'i 

^ EnUr any £. Urdtr DE. 

n 

" 

| 


DE : lx^2*y 1 1 +2x*y 1 -6y=G 1 


Exl : y 1 1 =x*y 

Indep. uar.=: lx 

Dep. oar.=: |y 

OPTIONAL: Initial condition 

at x or t on . . =: 12 

y or s or . . =: 110 

y 1 or s 1 or . . =: 115 

, CEnter=0K h <TESC=CFlNCEO ; 

MAIN RAD AUTD FUNC 1/30 


Interestingly this DE cannot be solved although it looks very similar to Example (4). Maybe that it 
doesn't have any solution? I try the 2. Order DE Checker offered in Option 4 foryi andy 2 : 


fl . Order 

F5T Yfhty FEt T FfiT 1 

lISsMaSOrder n|PDE|Transf |Moreh 1 

Differ 

1 : Any 2. Order DE 

2: Homogeneous 

3: Non- Homogeneous 

5! IUf So 1 uer : xl^x2 
6:Uariation of Para. 
7: Undetermined Coeff 
8: Bessel Equation 
9: Reduction of Order 
A: Legendre DEQ 
E: Cauchy-Euler DEQ 





MAIN RAD AUTD FUNC 1/30 
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’.-S-icWS *1 1 

i i 1 1 — 

-2-Cx-y' - 3 • y) 


y - 

X 2 


Thus, 

your solution: 


1 

y_ X 3 



is 

true 

Done. 



MAIN 

RAD AUTD 

FUhK 1/30 I7HH5 



y i and y 2 are solutions - as expected. So, we have solved the problem given in the textbook. I leave 
this DE for a later treatment and will look how DERIVE and other Computeralgebra Systems are per- 
forming. Doubtless I am starting again with DERIVE: 


# 2 : 


f 2 


D50LVE2_IV 


x 


2 

, 0, x, 2, 10, IS 

2 


■inappli' cable 


^ x J 

I was not very much surprised about #2. The TI’s algorithms are very close related with the DERIVE 
algorithms. (David Stoutemyer implemented the DE-package in DERIVE and the TIs as well.) 


I checked the validity ofyi andy 2 as solutions with DERIVE . Any linear combination ofyi andy 2 
should also form a solution, hence: 

2 -3 

sol (x) cl-x + c2-x 


2 

x ■ sol r r (x) + 2-x-sol r (x) - 6-sol(x) = 0 


So we are very curious what other systems will answer. Our next candidate is WIRIS: 

Example (5) 

prob5 : = solve(x 2 ■y ,, (x)+£x-y , (x)-6-y(x)=Q J y(2) = 10 J y , (2) = 15); 
prob5 — {Cy(x) = -16.e“ 3 - ln < K >+3.e 2 - ln < x >)} 

simplify(y(x) = -16e~ 3ln(x * +3e ?ln < x * ) — * y(x) = 3 X 16 

x 3 

We obtain for x > 0 y(x) = -16x -3 + 3x 2 . Then I try a compact form with wxMaxima. 

ic2 (ode2 (x A 2* T cliff (y, k, 2) +2*x* T cliff (y, x f 1) -6*y=0,y f x) , x=2,y=10, Miff (y f x) =15) ; 

? 16 

y = 3 x 

3 

x 

MuPad delivers the same solution without any problems. I don't know why DERIVE refuses solving 
this equation, do you know? 


I asked Albert Rich - one of the fathers of DERIVE - and I received his answer: 

Hello Josef, 

Thanks for your inquiry. Since David Stoutemyer wrote the DSOLVE packages for Derive, 
he would be the best one to help resolve the deficiency you found. 


Aloha, 

Albert 
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So my next mail was sent to David and he also answered very soon: 

Hello Josef! 

Karen and I are fine. We too took a trip to Tanzania this fall. Fantastic! 

DSOLVE2_IV(...) doesn't search for the pattern for Euler-type ODEs. 

However, often more than one method is applicable to a given ODE. One of the other imple- 
mented methods solved your first example, but not the second. 

However, it would clearly be a good idea to add Euler-type equations to the list of patterns to 
try. 

-- best regards, david 


Working out this review I remembered (supported by my old text books, of course) that DEs of this 
type - x 2 • y\x) + p • x • y\x ) + q • y(x) = r(x) - are Euler Equations - and fortunately enough, Nils has 
provided one option for solving this type of equations - even stepwise. DEQME in its recent version 
solves only the homogenous DE but we can use his tool for the inhomogeneous form, too. 

For the very few among you, who - like myself - have forgotten the algorithm solving an Euler DE, 
I’ll solve another example by hands and then check my calculations! 

Example (6) 

Find the general solution of x 2 y" - 3xy f + 8y = x 2 + 2x for x > 0. 

If y(x) is a solution for x > 0 then y{-x) is a solution for x < 0. 


The “trick” is applying the substitution x = e s . Then u(s) = y(e s ) =y(x) = u ( log x). 


The following is a useful application of the chain rule - for the students, 
x = e s and y(x) = y(e s ) = u(s) 


u'(s) = y'(e s )-e s =y'-x 


u\s ) = y"(e s )-e s -e s + y\e s )-e s =y"- e 2s + / . e* = / " . x 2 + / • x 
hence 

y\x)-x 2 = u\s)-u\s) and y\x)-x = u\s) 

The given differential equation reads now (in variables s und u)\ 
u"-u'-3u' + 8u = e 2s +2e s 
u"-4u' + 8u = e 2s +2e s 

This is an nonhomogeneous ODE of 2 nd order with constant coefficients. At first we have to find the 
general solution of the respective homogeneous DE. We ask DEQME\ 


f F1T 
|1 . Order 

KjBf FjT YFHtY F£t Y FfiT 1 

iWlISffianOrder n|PDE|Transf |More|* | 


1 s Any 2. Order DE 





Differ 

3: Non-Homogeneous 
4: 2. Order DE Checker 
5: IUP So Iyer: xl*x2 
6:Uariation of Para. 
7: Undetermined Coeff 
S: Bessel Equation 
0: Reduction of Order 
fl: Legendre DEQ 
B: Cauchy-Euler DEQ 



MAIN 


RAD liUTD FUhK £1/30 



r Y : :y Y "x'- Y i " Y i ■■■* 

i”Lv^L;r.i:|-: i\ z |0t^-;- ! C=|C L :i-= =i = 



r End Drdtr Homogeneous DiffEin. 



d E y^dx E +b+dy/dx+c*y=0 

or y 1 1 +b*y 1 +c*y=0 

has characteristic equation: 

s E+b* s+c =0 

Soloing it yields 3 different 
solution cases to the Deqn. 



^CEnter=0K ) CESC=CflHCEO J 




MAIN 

RAD AUTD FUNC £1/30 
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As we are restricted to y and x in this option, we 
have to rewrite the solution as 

u(s ) = c x • e 2s cos(2s) + c 2 • e 2s sin(2 s). 

We need a particular solution to accomplish the 
solution for the nonhomogeneous equation. 

The method of undetermined coefficients seems 
to be appropriate. Again we ask DEQME for 
information and support. 





Now we have solved the (u,s )- DE which must be transformed back to the y(x)-function: 
x = e s s = log x 


(u,s): u(s) = c x • e 2s cos(2s) + c 2 • e 2s sin(2s) + 


2J_ 

5 


(y, x): 


r(.v) -c x -x 2 cos(2 log x) + c 2 ■ x 2 sin(2 log x) + 



I regret having to admit that DERIVE is unable to solve this equation, so we will try with wxMaxima : 
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ode 2 (x A 2* T dif f (y f K f 2 ) t dif f (y f h, 1) + 8*y=x A 2 + 2*x, y, x) ; 
f 2 ) 2 ( 2 ) 2 

\5 x + 8 x) sin( 2 log( x ) ) + \5 x + 8 x) cos ( 2 log ( x ) ) , 

y = + x (§£lsin(2 1og(x)) + %k2 cos(2 log( x ))) 

20 

I used the occasion to demonstrate Options 2 and 3 of the F2-menu. 

We can apply the DEQME built-in DE Checker (Option 4) or do it on our own in the Home Screen 
(see below) for checking the solution. 


tm 


T , 4 l 

has general solution: 

, _ x 2 + x ■ ( 3 ■ y 1 + 2) - S y 


■=■ - 

X 2 


Thus, 

your solution: 


y = kl 
is 

true 

■x 2 cos(2- ln(x)) + k2 ■ x 2 ■ sin(2 

■ ln(x)| 

MAIN 

RAD AUTD FUNC 3/30 

I mss 


r= J ■■ || iii;Y fm . m r Y i : y Y Y _Y ^ 

□ 

:l 

r £. Urdtr DE Solution Checker. 


DE: : |3y Vx-8y^x^2+<x^2+2x )/x A 2 1 


H0TE1: Must be in format y''= 
N0TE2: Variables are x and y. 
Exl : y 1 1 =x+y 
Ex2: y 1 1 =y 

Solution y=: lkl*x^2*oos<:21n{x} !>-H 

Exl: y=x'"'2+x 

, CEnter=0K } CESC=CRHCEL4 ; 




MAIN RAD AUTD FUNC H/30 


ifl^T FF? Y F]:t T FHt y ft Y fT? Y '| 

F fr— |nigebra|Calo|0ther|FrgmI0|ciean Up| [ 


Done 


1 nightly = S3 ■ x 2 ■ cos(2 ■ ln(x^ + G4 ■ x 2 ■ sin> 
S3 ■ x 2 ■ cos(2 ■ ln(x)) + S4 ■ x 2 ■ sin(2 ■ ln(x)) > 


yZ.A 

dx 


^r(y) - 3 ■ x ■ ~ j-( y) + S ■ y x (x + 2) 
* dx 1 


x A 2*cK^ J .x J .2>-3x*t£<i; J .x J .l>+8^l 

MAIN RAD AUTD FUNC fi/30 


As you might guess I was not really satisfied about DERIVE' s inability solving one or the other Euler 
Equation. So look at this: 

(7) x 2 y" - 3x y f + 8y = x 2 + 2x 

(8) x 2 y" + 2xy’-6y = 0; y( 2) = 1 0, y' (2) = 1 5 

(9) x 2 / + 2xy’ -6 y = 30x; y( 2) = 10, y(6) = 15 

(10) x 2 y" + 2xy - 6y = 30x; y( 2) = 10, y'(6 ) = 15 

(1 1) x 2 / + 2xy'-6y = 30x; /( 2) = 1 0, y' (6) = 1 5 


All of them cannot be solved using DS0LVE2, DS0LVE2_IV and DS0LVE2_BV respectively. 
My EULER-functions do a better job. 

Examples (7) - (11) 


f 


DS0LVE2 


v 


3 

x 


S 


2 

x 


2 

x + 2 -x 


2 

x ^ 


inapplicable 


but now 


r 


EULER2 


v 


3 

x 


S 


2 

x 


2 

x + 2 -x 


2 

x > 


2 

2 2 x 2 ■ x 

cl ■ x -C0S(2 ■ LN(x)) + c2-x ■ SIN(2 ■ LN(x)) + + 

4 5 
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EULER2_IV 


0, x, 2, 10, 15 


x 2 

v x 


2 16 
3-x 


Parameter #9 in EULER2_BV (= k) = 0 (= default) for yO=y(xQ) and yl = y(xl) 

= 1 for yO = y(x0) and yl = y'(xl) 

= 2 for yO = y'(xO) and yl = y'(xl) 


EULER2_BV 


EULER2_BV 


EULER2_BV 


one check: 


6 

30 

i 

2 

X 

X 

6 

30 

l 

2 

X 

X 

6 

30 

l 

2 

X 


, x, 2, 10, 6, 15 


145 -x 15 -x 1620 


SS 


2 3 

11- x 


, x, 2, 10, 6, 15, 1 


310 -x 15 -x 2 26 SO 


163 


2 3 

163- x 


, x, 2, 10, 6, 15, 2 


1805-x 15 -x 64 SO 


96 S 


2 3 

121- x 


1S05-X 15- x 64 SO 


ffx) 


96 S 


2 3 

121- x 


2 

x -f rr (x) + 2 - x - f r (x ) - 6 - f (x) = HBH 

We will ask DEQME to solve example (8) after having recognized that (8) is of Euler type on 
page 43 investigating Option B: Cauchy-Euler DEQ. 

wxMaxima confirms the solution of boundary value problem (9): 


bc2 (ode2 (x A 2* T diff (y, x, 2 ) + 2*x* T dif f (y, x) - 6*y=3 0*x, y , x) , x=2,y=10, x=6,y=15) ; 
2 

145 x 15 x 1620 


muPAD says that I am right with my solutions for (9), (10) and (11), thank you very much, indeed. 

" solve (ode ({x A 2*y T T (x) +2*x*y T (x) -6*y (x) =30*x , y ( 2 ) =10 , y ( 6 ) =15 } , y (x) ) ) 

/ 145 jr 2 _ 15-i 1620 \ 

I 88 2 n.^3 / 


solve (ode ( { x 2*y T T (x) +2*x*y T (x)-6*y(x)=30*x,y(2) =10 ,y T (6) =15 } , y (x) ) ) 

f 310 jc 2 15 j: . 22680 \ 
l 163 2 163 x 3 J 


solve (ode ( { x A 2*y T T (x) +2*x*y T (x) -6*y (x)=30*x , y T ( 2 ) =10 , y T ( 6) =15 } , y (x) ) ) 

f 1805 x 2 15 -x 6480 \ 

I 968 2 12 l.x 3 J 


p 36 


Working with Nils Hahnfeld's DEQME (2) 


D-N-L#80 


These are two of my “home made” EULER-functions for DERIVE: 

EULER2 Cp , q i r, x, cl, c2, ans) : = 

Pnog 

ans := D5QLVE2(p-x - 1, q-x A 2, SUBST( r-x A 2 , x, e A x) , x, cl, c2) 

If - STRING?(ans) 

LIM(ans, x, LN(x)) 

"■inappl! cable" 

EULER2_IV(p, q, r, x, xO, yO, vO, ans , v, s) := 

Pnog 

ans := D5QLVE2(p-x - 1, q-x A 2, 5UB5T(r-x A 2, x, e A x) , x, cl, c2) 

If - STRING?(ans) 

Pnog 

v := 3(ans , x) 

s := (50LUTIQNS([LIM(an5 , x, xO) = yO, LIM(v, x , xO) = vO] , [cl, c2]))^l 
SUBST(ans , [cl, c2], s) 

11 -inappl! cable" 

I used these functions to solve problems (9) to (1 1). 


I don't copy the Boundary Value function EULER2_BV here to save space. Its syntax is as follows: 

EULER2_BV(p, q, r 5 x, xl, yl, x2, y2, k) 

with k = 0 (by default): given are (xi, y \ = y(xi)) and (x 2 , y 2 = y(x 2 )) 
with k= 1 : given are (xi, y\ = y(xi)) and (x 2 , y 2 = y’(x 2 )) and 

with k = 2: given are (xi, y\ = y’(xi)) and (x 2 , y 2 = y\x 2 )). 


WIRIS cannot solve nonhomogeneous Euler DEs. 
This is function euler2 for the TI-Nspire : 


j 2 2 ' 

deSoivey — y'+—y= CW 

t x X 2 j 

2 

y=c1 x +c2 x 

euler2 3/4 

Define euler2(/? J ^r J x J y)= 
Func 

Local kl ) k2 ) so!_ ) eq 

eq:=y"+{p-l)y'+qy=r\x=e l 
so 1_ : =de So 1 ve ( e q } ty) 
so/_k=ln(x) 

EndFunc 

de Solve 

2 6 1 

y+—y — -y=o,w 

x * 2 j 

„ ~2 , [x-y-3-y) 

x 2 

euler2[ : 2,2,0,x,y) y=c16x 2 +c15x 

eu/er2[2,-6.0,x l y) 2 c19 

' y=c20x A + 

x 3 


Nspire has the same problems solving Euler DEs as the TI-92 PLUS and the Voyage 200 as well. 
Please read the respective communication with Albert Rich and David Stoutemyer at the end of this 
article. 


Shame on me, I was not able to produce the I VP and BVP-functions for TI-Nspire because I din't find 
an easy way to solve the simultaneous equations for the constants cn (@n on the Tl-handhelds can be 
reset by pressing 8:Clear Home.). 


The screen shot on the next page shows some examples worked with TI-Nspire applying euler2. 
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eulerh 3, S,x 2 +2-x,x l y) y=c21x 2 - cos(2 ■ ln(x)) 

+c22'X 2 'Sin(2'ln(x))+d 5 

20 

de So Ive y 

y+--y=0 and y(l)=3 andy(l)=l,xy 

y=5'X-2'X 2 

\ 

* x 2 ] 


deSolve(y' 

-3- xy'+8-y=x^+2- x, x jv) 

y"=x 2 +X’(3y+2)-S^ 

© but, using euler2: 

eulerJ^-Z, S l x 2 +2-x,x,y) y=e3' x 2 ■ cos(2 ■ ln(x 

)Kc4.x 2 . s i„(21„(x)) + x(5X+8) 
; " 20 

© solving Problem 10: 

e n le v2{ 2j - 6j 3 0 ■ x, x^y) 

y=c6 x 2 - 15 ' X + CS 



2 X 3 

/(x):=right 

_ 2 Ib'X c5\ 
y=c6x + 

Done 

2 x 3 l 


so lve|l 0=_/(2) and 1 5 = (/(x) I x= 6, { c5, c6 } 

c 22680 J , 310 
c5= and c6= 

163 163 

310 2 

>- x 

15 -x 22680 
+ 

310'X 2 15'X 22680 

163 

2 163'X 3 

-| ■ | 

’ 163 2 'l63x 3 

1 


DEQME supports initial value and boundary value problems as well in Option 5: 


Diffet 


_ F]:t Y FHtT FFt YFfi^fT 

■Order n POE Transf More ► 


1 : Any 2. Order DE 
2: Homogeneous 
3: Non-Romogeneous 
4: 2. Order DE Checker 


Variation of Para. 
7: Undetermined Coeff 
8: Bessel Equation 
9: Reduction of Order 
R: Legendre DEQ 
B: Cauchy-Euler DEQ 



It works pretty well for all 2 nd order ODEs which can be solved by the TTs DE-Solver. 


Proceeding investigating menu F2 I reach Option 6: Variation of Parameters. DEQME provides 
stepwise demonstrating this standard algorithm. I’ll show two examples: 

(12) y" - 4 y' + 4y = 2e 4x 

(1 3) / - 9y = sinh(2x); y(l) = 1, y'(l) = -0.5 




( ?'. mv Y i i ^ Y Y #i‘ Y fii'r 

|4 1 


^ Vocation of PaKarotttK "> 


Solve y 1 1 +q£x>*y 1 +r{x}*y=g{x> 
If y 1 1 +q(x)*y 1 +r t x>*y=0 has 
solutions yl(x) and y2tx> 
then partic. solution is: 
VCx>=vl(x>*yKx>+v2Cx>*y2(:x> 
where vl(x>= 
-Yy2*gx{yl*y2'-y2*yl 1 >dx 
v2=Y(yl*gx<yl*y2'-y2*yl 1 >dx 

, 6Enter=0K } 6ESC=CRNCEL} ; 





MAIN 

RAD AUTD FUNC 0/30 





Maybe that students are happy with this, notify the result and proceed to the next example. I am not 
happy, because I can see that both integrals are not calculated. So this method fails! At the other hand 
there is a result, which says that the internal DE Solver finds the particular solution! What is my con- 
clusion? What should the students conclude? There must be another method which should be more 
appropriate. Try in this case the Method of Undetermined Constants! 


DERIVE returns - without stepwise simplification: 

4-x 

4-x e 2-x 

D5QLVE2C-4, 4, 2-e ) = + s ■ (c2-x + cl) 

2 


I recommend applying DERIVE' s Stepwise Simplification of this Solving procedure, it is an experi- 
ence following the 30 “steps”. You will learn a lot about DERIVE ' s “thinking”. 


I am trying Parameter Variation once more solving example (13): 
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HII'Y Y : • r Y "w Y f" Y t-.-v 

n 


r Variation of Pararotttr 


DE<x,y>: |y "-9y=sinh<2x> 1 


Exl : y ,l +y=tantx> 

Coeff of y 1 1 nust be troade} 1. 
OPTIOHRL: Initial condition 

at x=: lln<2> 

y=: 11 

y'=: 1-1x2 

, CEnter=0K <TESC=CFlNCEO, 




MAIN RAD AUTU FUNC 0/30 



:i'!!!!: 1^ L L3 -.eJ^: : E= : =1: ^ j C k : 


-te-'i -5*x5*Cp''C4*x)-1 >*C^(6*x> + l V4x(6> 
-J r<? A < -2*x>*<e'-<2*x> + l >*(e A (6*x> + l V<12*f 

e’ 2 *-^ -J3 ■* 2 - x -tan-i('^ 2 y 2 * + ^ 

■)-> 

72 


and u2(x>= 

Tyl*g/(yl*y2'-y2+yl 1 >dx 


MAIN RAD AUTD FUNC 0/30 

17HH3 


fi/ .j!!!!!: [}-u : ^T3-.=:1 

m 




□ 

2> Particular solution 

ul(x>*yl(x>+o2(:x>+y2t:x> = 

u2(x) e " 3x (e 6x + 4e 3x + l) 

nh(3*x)-co^ 

x)+l?x4x(6f 

2 

where vl(x)= 

-Ty2*qX(;yl*y2 1 -y2*yl 1 )dx 
"Tcosh < 3*x ) *s i nh < 2*x > ^ < 2*3*s i 
-re-'- < -5*x>*<* A <:4*x>-l >*<* A <6* 

MAIN RAD AUTD 


FUNC 0/30 IdVUHa 



and o2(x>= 

J-yl*q/{yl*y2'-y2*yl 1 >dx 

T2*s i nh < 2*x V 2*3*s i nh < 3*x > -cosh 3*x > *0^ 

T 2*s i nh < 2*x !> / < 6*s i nh < 3*x > -0 > dx= 

^ A x+Cf A (2*x)+l)/(3+Cf A (4*x)+f A C2*x)+l)^ 



9 

MlilN Fill [i liUTD FUNC 0/50 IMIfeH 



Thus, Particular Solution 

ul(x>*ul(x>+u2(:x>+u2t:x> = 

e ~ 3 x (l0-e4f 6 x - g 5 * + g* + 10 63) 


10 

Using Initial condition, 
particular solution is: 

e ' 3 x (l79 f 6 j< - 192-e 5 '* + 192 e x + ■. 
t ~ 1920 

MlilN Fill [i liUTD FUNC 0/50 imfeH 


Compare with DERIVE followed by wxMaxima : 

D50LVE2_IV(Q, -9, SINH(2-x), x, LN(2), 1, -0.S) 


3x 

179 -e 


1920 


2-x - 2-x - 3-x 

e e 151- e 

+ + 

10 10 30 


ic2 (ode 2 (Miff (y f x, 2 ) - 9*y=sinh (2*x) f y f x) f x=log (2 ) , y=l f Miff (y f x) =- 1/ 2 ) ; 

- 2 x { 4 x H '| 3 x - 3 x 

“1/ 179 % e 151 %e 

y = + + 

1G 1920 30 

F2 Option 7: Undetermined Coeff(icients) was treated on page 33. 

Options 6 and 7 are really nice and students and teachers as well might appreciate this tool. It could 
be a challenge to implement this “stepwise simplification” for TI-Nspire or DERIVE or other Comput- 
eralgebra systems - as stand alone programs. 

So we can proceed to Option 8: Bessel Equation. 

DEQME gives a short explanation of how the Bessel Equation looks like. The general form of the 
Bessel DE is: 

x 2 y" + xy f + (x 2 -m 2 )y = 0. 

I learned from the text books that the Bessel DE and Legendre Equation (Option A) appear when 
solving Partial Differential Equations ( Laplace Equation). 

Obviously DEQME shows only the - most important - case for m = 0. 
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Bt55fcl E -motion. 


x E *d E yxdx 2 +x*dyxdx+x 2 *y=0 

or the series solution: 

y=C*2<< -l^n*x"(:2n>>x(;2-"<:2n>*(:n! > s > 

■:: ;Enter=0K > CESC=CflHCEL: > 


MAIH FiliD AUTD FUhK E/30 


SeeMuPad’s answer: 

solve (ode (x 2*y T T (x) +x*y T (x) +x A 2*y (x) =0 ,y(x))) 

Warning: Only Q-solvable exponential solutions will be found! [ode : secondOrder] 
(C2-J 0 (x) + C3 -¥„(*)} 

and wxMaxima ' s answer: 


(%i3 0) ode2 (x A 2* T dif f (y, x, 2) +x* T diff (y, x) +x A 2*y=0, y, x) ; 

( %o3 0 ) y = bessel_y ( 0 f x ) %k2 + bessel_j ( 0 f x ) %kl 

(%i31) ode2 (x A 2* T dif f (y, x, 2) +x* T diff (y, x) + (x A 2- 4 ) *y=0 , y , x) ; 

(%o31) y = bessel_y ( 2 f x ) %k2 + bessel_j ( 2 f x)%kl 

( % i 3 2 ) ode2 (x A 2* T diff (y, x, 2) +x* T diff (y, x) + (x A 2-3 ) *y=0 , y , x) ; 
is an integer? type y or n. 

The solution contains the so called Bessel functions (see DNLs #18 and #34). DERIVE provides a 
utility file BesselFunctions .mth . Bessel functions can also be found in SpecialFunc- 
tions . dfw and SpecialFunctions . mth provided in the Users Folder. 


After this short trip into the world of higher mathematics I will return to easier fields. 


When I found F2 Option 9: Reduction of Order I - again - consulted my text books on ODEs and 
didn't find anything matching with this option. I didn't find one single problem similar to (14) or (15). 
Given is a 2 nd order ODE and one of its solutions. Find the second one. 

(14) y”(l-x 2 )-2xy' + 6y = 0; y x =ff 

( 15 ) y" - 2y r - 3y = 6; y t =e 3 * 


fl . Order 

F]:t YFHtY FEt Y FbT ff/ } 

311ffia30rder n|PDE|Transf |Moreh | 

Differ 

1 : Any 2. Order DE 
2: Homogeneous 
3: Non-Romogeneous 
4: 2. Order DE Checker 
5: IUP Soluer: xl*x2 
6:Uariation of Para. 
7: Undetermined Coeff 
8: Bessel Equation 


R: Legendre DEO 
B: Cauchy-Euler DEQ 




MAIN RAD AUTD FUhK £3/30 
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2. Solution y2<x>= 


yl<x>*,Me A t -XqCx>dxVyl(x>2 ]dx 
Here, - Xq{x>dx= 

S -2x^Cl-x A 2)dx= 

J'-2*xx(x- A -2-l>dx= 

■ 1 n( x 2 - l) 

and C -J'q(:x>dx> = 

e--:: -lnCx"2-l>dx> = 

MlilN Fill [i liUTD FUNC £5/50 \:m>\ 


rjjskvsu 

■ ■ v "l 

i I c | 

o&v 


fc,,;s *1 

□ 

r i r id 1 1 y , — yzA J 

4l(x)*Uf A ( ", 
(3x^2-l >x2*J' i 
(3x^2-lV2+-: 

-(3 x 2 - l) ■ li 

!T7= 

rciCx)dxVul<x)2 ]dx= 

<x^2-l >*{3*x^2-l > A 2>dx= 
InCx+l >x2+lnt!x-lV2+r C3>/(:2^ 

^x +1) . (3 ■ x 2 - l)- ln(x - 1) 

4 

(l - 3-x 2 )- InCx + 1) +(3-; 

< 2 -0- 

4 

InCx - 1) + ^ 

4 

MAIN RAD AUTU 


FUNC £3/34 IdVUHa 



I check the validity of the complete solution - and I can be satisfied. DEQME does a good job. 


3-x - 1 

loes(x) := cl + c2- 


(3-x - 1) ■ LN 


x - 1 


x + 1 


3-x 


2 J 


loss' p (x)-(l - x ) - 2-x-loes p (x) + 6-loes(x) = 0 


Example (15) 




^CEnter=0K > CESC=CflHCEL^ 



E :4 ) 

e 2 * 

Thus, 

-^(x)dx)/Ml(x)S]dx= 

Xf A C2*x)/e A (3x) A 2= 


- e -4x 


4 


Finally, y2{x!> = 

yl(x>*T[^C -TyCx>dxVyl(x>2 ]dx= 


MAIN RAD AUTD FUNC 1/34 

liiTIIHa 


The complete solution of the homogeneous equation is c x e 3x + c 2 e x . Check it! 

I asked Nils about the relevance of this option because I couldn't find any related problem in my 
books. He answered: 

Yes, “Reduction of order” is described here 
http://tutorial.math.lamar.edu/Classes/DE/ReductionofOrder.aspx 


Is treated in all good ODE courses. 

(This is really a fine web site. Try this URL, Josef) 
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I mentioned the Legendre Equation earlier. Option A treats this kind of 2 nd order ODE. 
(16) Legendre DE: (l-x 2 )y" -2xy' + «(« + l) = 0; neN. 


fl . Order 

F5T Y F4t Y FEt Y FfiT 1 

lISsMaSOrder n|PDE|Transf |Moreh 1 

Differ 

1 : Any 2. Order DE 
2: Homogeneous 
3: Non-flomogeneous 
4: 2. Order DE Checker 
5: IUP Soloer: xl*x2 
6:Uariation of Para. 
7: Undetermined Coeff 
8: Bessel Equation 
9: Reduction of Order 




fi: Cauchy-Kler 







MAIN 

RAD AUTD FUNC 1/30 




ffiT^OY Flnr Y f::t Y FHt Y fe Y FfiT 

n 


f DiFH E-in. "‘i 

t l-x A 2>y 1 l -2xy , +n*(n+l)y=0 

n: 121 | 

Choose l^n^5 . 

Solution is a*yl(x)+b i t ( y2(x) with 
yl tx)=Legendre Polynomial. 
y2£x) results from Reduction 
of Order. 

k CEnier=0K ) CESC=CFINCED, 

IwO i 

MAIN RAD AUTD FUNC 0/30 




1 

(3 ■ X 

2 - l)- ln(x +1) [3 x 2 - l)- 

ln(x-l) h 


4 4 


-(Ci 

- 3x 2 )- InCx + 1) +(3 x 2 - 1 

llnCx- 1)> 

Thus 

b (3 

4 

, a*yl tx>+b*y2tx>= 

■x 2 - l)- InCx + 1) b (3 ■ x 2 

- 1 ) • 1 n( x - ► 


4 

4 
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Please compare with (14) from above. Now you might find an explanation for my choice of v, . 
The Legendre Polynomials for n e N can be produced by the “Formula of Rodrigues”: 

l fd y ’ 


L in) =- 


2" n\ 


dx 


(* - 1 )". 


VECTOR 

' 1 

'd ^ 

k 2 k ^ 

(x - 1) , k, 0, 5 

k 

5.dx^ 




L 2 -k! 


J 


2 2 4 2 4 2 

3-x -1 x-(5-x -3) 35-x - 30-x +3 x-(63-x - 70-x +15) 

1 1 i i i 

2 2 8 8 

As calculation times for finding the Legendre polynomials increase enourmously on the TIs, Nils re- 
stricted for n < 5. 


MuPAD : 

solve (ode ( (l-x A 2) *y T 1 (x) -2*x*y 1 (x) +2*y (x) =0, y (x) ) ) 

{ C6 .* + cs - A22£^ilL + 1) } 

solve (ode ( (l-x A 2) *y T T (x) -2*x*y T (x) +6*y (x) =0,y (x) ) ) 

|c3- + C2- ( — 3 h t ;f " 1 ) | 3 ln t :f+1 ) , 9-jf 3 -h(jr+ 1) 

( \ 3 / \ 4 3 3 3 3 

(Nice job for the students: Compare the result with the result from page 40, Problem (14)!) 


By the way, the LDE - algorithm works also for n = 0. 
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The last option doesn’t need much space, because we talked about 2 nd order Euler Equations and how 
to treat them on page 32. Option B does this in a very reduced form See example (8) treated with this 
option: 


fl . Order 

F3t YFHtY Y FbT ’p -1 

IffllfllgHaQrder n|PDE|Transf |More|* | 

Differ 

1 : Any 2. Order DE 
2: Homogeneous 
3: Non-flomogeneous 
4: 2. Order DE Checker 
5: IUP Soluer: xl*x2 
6:Uariation of Para. 
7: Undetermined Coeff 
8: Bessel Equation 
9: Reduction of Order 
fl: Legendre DEQ 


MAIN RAD AUTD FUNC 1/30 


ESP 


£. Order Couchy-Eulth- DEC. 


„nl 

Euler DEO. > ^-1 > 


haoe the format.: 
x^2*y 1 1 +b*x*y 1 +c*y=0 
Solution: 

using x=e^s} we get: 
y 1 '+tb-l>*y '+c*y=0 
which requires us to soloe 
s^2+tb-l )*s+c=0 


'■| Enter=OK > <ESC=CnRCEL> 



f 3 solution costs. 

1 

O If <b-l>2-4c>0, then 
y—C 1 < s 1 *s > +C2*e^ < s2*s 
2> If <b-l>2-4c=0 ? then 
y=[Cl+C2*ln<x> ]*x A (sl) 

3!> If <b-l>2-4c<0, then 
y=[Cl*sin<E*lntx>> 
+C2*cos<E*lntx> ') ]*lxl A fl 
where fl= _ { b- 1 ;■ .-'2 , 
E=lx2+fC4c-Cb-n^. 

, CEnter=0K } tESCCflNCEO ; 
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I am looking forward to investigating Menu F3 
Order n containing the topics Linar+Const 
Coeff, 2xLinear System, X’ = A * X, 

X’ = A * X + F and Separable DE - and to 
“remember” a lot about Differential Equations 
again. 


Differential Equations Made Easy for TI-89, 92+ and Voyage 200 is available from 

http://www.ti89.com 


You can find much information about Legendre and Bessel DEs at 

http://mo.mathematik.uni-stuttqart.de/inhalt/beispiel/beispiel822/ 

http://en.wikiversitv.orq/wiki/Leqendre differential equation 

http://www.math.tuqraz.at/~berqlez/Math C/Folien 1 0%28Potenzreihenans%29.pdf 
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Differential Equations “stepwise” with DERIVE and Tl-Nspire 

Josef Bohm, Wiirmla, Austria 

I liked the stepwise simplification in DEQME - and Nils wrote that the users of DEQME ap- 
preciate this feature of DEQME, too. On page 39 I mentioned the “challenge” how to create a 
similar procedure with DERIVE or Tl-Nspire or any other CAS-tool. 

Nils and I, we had an extended exchange of emails concerning some features and improve- 
ments of the F2 Menu options. So I found some time to face my own challenge. Here are my 
results for the homogenous 2 nd order ODEs with constant coefficients: 

The DERIVE Stepwise (without printing the progam): 


"Stepwise solution for homogeneous 2nd order ODEs with constant coefficients, 
y" + b y 1 + c y = 0 

hom2JV(b 1 c 1 x0 1 y0,v0 1 x,y) with y(xO)=yO and y'(xO)=vO. 

Three examples: 

y" - 4y = 0 

#1: horn2_IV(0, -4) 

Characteristic equation: 
s = -2 v s = 2 
Discriminant = 16 

Discriminant > 0 ■+ two real solutions: si £ s2 

[2 i -2] 

General solution: y = cl*e A (sl*x)+c2*e A (s2*x) , 
hence : 

2- x - 2x 

#2 : y = cl-e + c2-e 

x M (t) - 6 x'(t) + 9 x(t) = 0 

#3: hom2_IV(-6 , 9, t, x) 

Characteristic equation: 
s A 2 - 6-s + 9 = 0 
Discriminant = 0 

Discriminant = 0 -+ one real double solution: si = s2 

[3] 

hence 

3- t 3-t 

#4: x = cl-e + c2-t-e 
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y" _ 4y' + 8y = 0, y(7T/2) = 1 , yV2) = -1 


#5: 


hora2_IV 



s, X, y, 




Characteristic equation: 
s a 2 - 4 ■ s + S = 0 
Discriminant = -16 

Discriminant < 0 two complex solutions: si ± s2 
[2 + 2 - l , 2 - 2 ■ i_] 

General solution: y = cl*e A (RE(sl)*x)wsin(IM(sl)*x)+c2*e A (RE(s2)*x)*cos(IM(s2)wx) , 
hence : 

General solution: 

y = e A (2-xMc2-C0S(2«x) + cl-SIN(2-x)) 

Solving the Initial Value Problem: 

1. equation: y(x0)=y0 
c2 = - e A (-7r) 

2. equation: y(x0)=v0 
2 -cl + 2 ■ c2 = e A (-tt) 

Solutions for cl and c2 are: 

[[3.e A C-^)/2], [- e A (-jO]] 

Solution 


# 6 : 


# 7 : 


2-x 

General solution y = e ■ (c2-C0S(2-x) + cl-SIN(2-x)) 

2-x - n f 3- SIN(2-x) 

Therefore, the solution is: y = a ■ - C0S(2-x) 

l 2 



^ n ^ 

2-X - TT 

f 3- SIN (2-x) *1 

DS0LVE2_IV 

-4, S, 0, x, ,1,-1 

= e 

C0S(2-x) 


L 2 J 


l 2 J 


The “steps” printed in blue are created by using the DISPLAY-function. We don’t have a 
PAUSE-function like in the TI-89, 92 and Voyage 200 programming language. So I can show 
the steps in a form of report only. 

Unfortunately I cannot enter “DERIVE s interior” in order to make use of DERIVE ' s stepwise 
calculation. As I mentioned ealier it is very interesting to apply the DERIVE steps on the 
DERIVE functions. Sometimes you might fail. I entered DSOLVE2_IV(-4, 8, 0,x,tc/2, 1 ,-1 ) and 
was stopped by “Memory exhausted” after three steps .... 

I exchanged some mails with Jose Luis Galan from Spain. He wrote that he and his col- 
leagues are working with this type of functions (including explanatory comments) on their 
university. 


p 46 


Josef Bohm: 2 nd order DEs “stepwise” 


D-N-L#80 


Using the features of the latest version of Tl-Nspire (there are rumors that release of version 
3 will be soon) results also in a nice dialogue-driven commented output of the algorithm. 

The Request-, Text-, and Disp-command are of importance. The values for b and c, and for 
xO, y(xO) and y’(xO) are entered in a dialogue box. 
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2. Equation: y'(xG)=vQ 


Solutions for cl and c2 are: 


■ and c2=-^~ 


cl- 

2-e n e 71 

Special solution is 


3'(e x )" J 'sin(2'x) (e x ) ■cos(2 , x) 


2-e 71 


Done 


The screen shots show the Calculator Application. 


This is the first part of the Tl-Nspire program. 

Below is the start of the DERIVE function. 

Both programs can be downloaded from the 
DUG-website (contained in mth80.zip). 


hom2_iv 

Define hom2_iv(|)= 

Prgm 

Loca 1 c i vs } d, so Is } ge n_so I t y_ l ) eql i eq2 ) cs } xO,yO, vO 
Disp 11 DE of form y" + b*y‘ + c*y = 0" 

RequestStr 1 ' Enter b, c : " f c o 
co: = "{"&co&" }" :co:=expr(co) 
fi:=co[l]: c:=co[2] 

Disp "hom. DE: " 

Disp y"+by'+cy=Q 

Disp "Enter In itia lvalues xO, y(xQ), v0=y'(x0)" 
Disp "ENTER 0 if none" 

RequestStr "In itia lvalues x0,y(x0),y'(x0):" J ivs 

?vs: = "{ M &?vs&" }" 

ivs:=expr(ivs) 

d\=b 2 -A-c 

so Is : =cZeros (s 2 + b • s+ c, s) 

Disp "Characteristic Equation: 11 

Disp s^+b-s+c=Q 


hom2_IV(b, c, x, y, xO := i. , yO, vO, d, s_, sols, gen_sol , sp_sol , cl, c2, 
Pnog 

d := t^2 - 4-c 

sols := S0LUTI0NS(s A 2 + b-s + c = 0, s) 

DISPLAY( M Character!sti'c equation : 11 ) 
s_ := P (s A 2 + b-s + c = 0) 

DISP LAY CO 

"DISPLAY(s A 2 + b-s + c = 0) M 
DISPLAY( M Discriminant" = d) 

If d > 0 
Pnog 

DISPLAY(" Discriminant > 0 -* tijo real solutions: si t 52") 
DISPLAY(sols) 

DISPLAY("General solution: y = cl*e A (sl*x)+c2*e A (s2*x) , ") 
DISPLAYC'hence: 11 ) 

gen_sol := y = cl-e A (x- sols.j.1) + c2-e A (x- sols.j.2) 

If d = 0 
Pnog 

DISPLAY( M Discriminant = 0 -* one real double solution: si = s2") 
DISPLAY(sols) 

DISP LAY ("General solution: y = cl*e A (sl*x)+c2*x*e A (s2*x) , ") 
DISPLAYC'hence: 11 ) 

gen_sol := y = cl-e A (x- sols.j.1) + c2-x-e A (x- solsRl) 
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Dear Noor and Josef, 

Great to hear from you! From the weather news I have been hearing from Europe, I can as- 
sume the picture of your village below is not a current one. :=) 

As far as my recent work is concerned, my goal is to pass on the lessons learned in my ca- 
reer to the next generation of CAS developers. To accomplish that, my website promotes 
using a rule-based approach to implement such systems. As proof-of-concept, the site 
makes freely the rules required to integrate a large class of expressions, and provides test 
results favorably comparing Rubi, a rule-based integrator, with the major commercial sys- 
tems. 

I am heartened to hear the Derive User Group continues to be the virtual home for the com- 
munity of loyal Derive users around the world. I only hope some entrepreneur will use the 
knowledge on my website to produce a worthy successor to Derive for the future genera- 
tions... 

Aloha, Albert 


A DUG Member from Germany has problems with the program editor window. You can see the Ok 
and the Abbrechen (= Cancel) button within the written text. Is there anybody facing similar prob- 
lems? Do you know how to get rid of these nasty buttons? 


RETURN DELETE (vNluin, nPos) 
vNum 


VekDelObjSmMax(xMax, vNutn, nPos := oo) : 
If nPos = co 
Prog 

nPos := DIM(vNum) 

Loop 

If nPos = 0 
RETURN vNum 
If vNurrunPos < xMax 

vNum := DELETE (vNum, nPos) 
nPos 1 

Prog 

If nPos > DIM(vNum) 

RETURN vNum 
If vNurrunPos > xMax 

RETURN DELETE CvNum, nPos) 
vNum 


Funktionsdefinition bearbeiten 

|x| 

Funktion: Name und Argumente: 

|VekDe 10b jSmMax (xMax , vNum, nPos := oo) 

3 

Funktionsdefinition: 



IF(nPos = oo, 

PR0G( 

nPos := DIM(vNum), 

L00P( 

IF(nPos = 0, 

RETURN vNum), 

IF (vNurrunPos < xMax, 

vNum := DELETE (vNum, nPos)), 
nPos 1) 

PROGC 
IF(nPos > DIM (vNum) , 

RETURN vNum) 

IF(vNununPos > x OK | Abbrechen 
RETURN DELETL'Cv^um~riPos) , ' 


), 


Do you know this? It is possible to include internet links in a DERIVE-file. Simply write 
the URL in a text box. It will appear blue and underlined and you can open the URL from 
within DERIVE. Josef 

( jt 2-x - n ( 3-5IN(2-x) "j 

#1: DSGLVE2_IVI- -4, 8, 0, x, — , 1, -ll = e -I C0S(2-x) j 

This is a comment within a textbox. 

You can immediately open the DUG website by clicking on 

www.austromath.at/duq 

or to any other website 

www.ti-nspire.com 

2-x 

#2: DSOLVE2C-4, 8, 0) = e ■ (cl-C0S(2-x) + c2-SIN(2-x)) 


Hello Derivers, the routines provided in the LinearAlgebra.mth file leave something to be desired in 
practice - has anyone resolutely attacked the problem of diagonalization and determination of eigen- 
spaces of large matrices? For my application large = 8x8. 1 need a robust and fast diagonalization 
scheme. 


Danny Ross Lunsford [antimatter33@YAHOO.COM] 




